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ABSTRACT 

The cosmic-ray positron and negatron spectra between 

11 and 204 MeV have been measured m a series of 3 high- 

altitude balloon flights launched from Fort Churchill, 

Manitoba, on July 16, July 21, and July 29, 1968. The 

detector system consisted of a magnetic spectrometer 

utilizing a 1000-gauss permanent magnet, scintillation 

v 

counters, and a lucite Cerenkov counter. 

Launches were timed so that the ascent through the 
100 g/cm 2 level of residual atmosphere occurred after the 
evening geomagnetic cutoff transition. Data gathered 
during ascent are used to correct for the contribution of 
atmospheric secondary electrons to the flux measured at 
float altitude. All flights floated near 2.4 g/cm 2 
residual atmosphere throughout the nighttime interval. 

A pronounced morning intensity increase was observed 
in each flight. We present daytime positron and negatron 
data which support the interpretation of the diurnal flux 
variation as a change in the local geomagnetic cutoff. A 
large diurnal variation was observed in the count rate of 
positrons and negatrons with magnetic rigidities less than 
11 MV and is evidence that the nighttime cutoff was well 
below this value. 
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Using nighttime data we derive extraterrestrial 
positron and negatron spectra. The positron- to- total- 
electron ratio which we measure indicates that the inter- 
stellar secondary, or collision, source contributes 
£50 percent of the electron flux within this energy 
interval. By comparing our measured positron spectrum 
with the positron spectrum calculated for the collision 
source we derive the absolute solar modulation for posi- 
trons in 1968. Assuming negligible energy loss during 
modulation, we derive the total interstellar electron 
spectrum as well as the spectrum of directly accelerated, 
or primary, electrons. We examine the effect of adiabatic 
deceleration and find that many of the conclusions 
regarding the interstellar electron spectrum are not 
significantly altered for an assumed energy loss of up to 
50 percent of the original energy. 
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I. 


INTRODUCTION 


Measurements of the composition and energy spectrum 
of cosmic rays are of importance to the study of a variety 
of astrophysical problems. The total energy content of the 
cosmic rays is sufficient to give them an important role in 
the dynamics of the galaxy. Knowledge of the galactic 
spectra of the cosmic-ray constituents is necessary to any 
theory which attempts to describe that role, as well as to 
theories of the origin, propagation, and containment of 
these energetic particles. The cosmic-ray spectra measured 
near the earth differ considerably from the interstellar 
spectra, however, due to the modulation effects of the 
solar wind. Measurements of the positron and negatron* 
spectra can provide a valuable tool for the study of the 
solar modulation of the cosmic-ray spectra. The only sig- 
nificant source of cosmic-ray positrons with energies above 
a few MeV is believed to be the decay of pions produced in 
collisions of high-energy cosmic-ray nuclei with the inter- 
stellar matter. The total solar modulation of the positron 
spectrum can be determined therefore by comparing measure- 
ments made near the earth with calculations of the 


*In this thesis the designations "positron" and 
"negatron" will be used whenever the sign of the electron 
charge is relevant to the discussion. 
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interstellar spectrum of secondary positrons. The measured 
electron charge ratio provides information on the relative 
contributions of electrons resulting from pion decay and 
from other sources, and, when combined with the calculated 
secondary positron spectrum, enables us to reconstruct the 
interstellar electron spectrum from all sources. Coupled 
with surveys of the non-thermal galactic radio emission, 
knowledge of the interstellar electron spectrum can also 
yield information about conditions in interstellar space, 
m particular about the magnetic field and the location and 
extent of H^-j regions. 

Electrons represent only a small fraction of the 
cosmic rays incident on the earth. Early attempts to deter- 
mine the primary electron flux (1-3) succeeded only in 
setting upper limits. The experiment of Cntchfield, Ney, 
and Oleksa (3), for example, set an upper limit of 0.6 
percent for the ratio of the flux of electrons with energy 
above 1 GeV to the total measured cosmic-ray flux. The 
presence of cosmic-ray electrons in the galaxy was deduced, 
however, from the observation of the galactic non-thermal 
radio noise. In 1950, Kiepenheuer (4) suggested that the 
galactic radio emission was synchroton radiation from high- 
energy electrons spiralling in the galactic magnetic field. 
He estimated that an electron flux of only one percent of 
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the cosmic-ray proton flux would be sufficient to account 
for the observed radio intensity. In 1952, Hayakawa (5) 
fxrst examined the electron flux resulting from the decay 
of charged pions produced in interstellar nuclear inter- 
actions. He found by a rough calculation that a flux of 
the magnitude required by Kiepenheuer could be entirely 
accounted for by the pion-decay source alone. 

The first direct observations of an extraterrestrial 
electron flux were made in 1960 by Meyer and Vogt (6) and 
by Earl (7) . The cosmic-ray electron spectrum at the 
earth has been investigated subsequently by many workers 
(8-21). (These references are primarily recent results. 
References to earlier publications will be found in many of 
these papers . ) 

Several authors, taking advantage of improved know- 
ledge of the relevant nuclear and galactic parameters, have 
recalculated the galactic equilibrium electron spectrum re- 
sulting from pion decay (22-29). The spectrum of interstel- 
lar knock-on electrons has also been calculated (30) and is 
found to dominate the pion-decay source below about 20 MeV. 

l 

These two sources together constitute the collision, or so- 
called "secondary", source of interstellar electrons. 
Alternative sources which have been suggested for the 
origin of cosmic-ray electrons (such as Fermi acceleration 
or acceleration in super-novae or pulsars (25) (31)) all 
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involve direct acceleration and hence are referred to as 
"primary" sources. The secondary source contains a high 
proportion of positrons at all except the lowest energies. 
The primary source, whatever its exact nature, very pro- 
bably injects only negatrons, since the energetic particles 
are accelerated out of the ambient matter, which presumably 
contains few free positrons. Hence a comparison of the 
electron charge ratio measured at the earth with that 
calculated for the secondary source alone permits the 
determination of the relative contributions of the primary 
and secondary cosmic-ray electron sources. 

Experiments to measure the electron charge ratio 
have been performed by several investigators. A University 
of Chicago group flew a balloon-borne magnetic spectrometer 
during the period from 1963 to 1966 (32-35) . The instru- 
ment covered an energy range from about 200 MeV to 14 BeV. 
Their measured positron spectrum agrees with calculations 
for positrons from pion decay but indicates that the major- 
ity of particles of these energies are of primary, rather 
than secondary, origin. Two other balloon-borne experi- 
ments have utilized the east-west asymmetry of the 
geomagnetic cutoff to determine electron charge ratios. A 
Saclay-University of Milan group flew a spark chamber ex- 
periment at a location with a vertical cutoff rigidity of 
5.4 GV (36) (15). They reported a positron fraction, 
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N 4 ./ CN , + N _) , of < 0.3 9 ± 0.11 between 4.6 and 6.5 GeV. 
e +/ e 

Daniel and Stephens flew oriented emulsions where the ver- 
tical cutoff rigidity was 16.2 GV (37). They claimed to 
see a positron fraction of 0.70 ± 0.20 between 15 and 50 
GeV. In later publications the quoted error limits have 
been considerably increased, however, and analysis of addi- 
tional events has indicated a possible negative excess (38) 
(39) . Their reported electron spectrum is considerably 
above the calculated secondary spectrum at these energies. 

Cline and Hones, using a satellite-borne detector, 
have reported a flux of 200 positrons/ (m 2 ♦ sec • sr* MeV) 
between 0.5 and 3.0 MeV (40). Their instrument identified 
positrons by detecting the annihilation photons and did not 
measure the total electron flux. Because of instrumental 
uncertainties the measurement was quoted as an upper limit. 
This very high positron flux is several orders of magnitude 
above the pion-decay source at these energies. Beta decay 
of radioactive spallation products in the cosmic rays has 
been examined as a possible source for these positrons (41) 
(42) . This source requires a very large interstellar flux 
of cosmic-ray nuclei, however. 

Recently Cline and Porreca (43) and Kniffen, Cline, 
and Fichtel (44) have reported measurements of the positron 
and negatron spectra between 2 and 10 MeV and between 35 
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and 220 MeV, respectively. We shall discuss their measure- 
ments m Chapter V. 

In this thesis we report measurements of the positron 

and negatron spectra between 11 and 204 MeV made during 

three high— altitude balloon flights launched from Fort 

Churchill, Manitoba, in July, 1968. We have concentrated 

on this energy interval since it was previously only poorly 

explored due to inherent experimental difficulties and was 

potentially of great value for solar modulation studies. 

Other means of studying the absolute solar modulation, 

such as measurements of the deuterium and helium-3 spectra 

or inferences from the galactic synchrotron radiation, are 

confined by experimental limitations to rigidities above 

several hundred MV. Previous measurements of the electron 

/ 

spectrum between 10 and 200 MeV at the earth fell well 
below the calculated interstellar secondary spectrum alone. 
Over most of this energy interval the calculated positron 
fracti'on of the secondary source is considerably greater 
than one-half. Since the data of the University of Chicago 
group indicated that the measured positron fraction might 
be increasing below ■ — 400 MeV, it appeared that a measur- 
able positron flux might be found below 200 MeV. Such 
measurements would yield valuable new information about 
solar modulation and the interstellar electron spectrum. 
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Measured electron spectra which overlap all or part 
of our energy range have been published previously (11) (18) 

(21) (45) (46). Although this reference list is confined 

to satellite experiments and to balloon experiments for 
which the local geomagnetic cutoff is known to fall below 
the detector threshold, there are nevertheless large dif- 

■» 

ferences m the reported fluxes. The variations may be 
attributable in part to solar modulation effects. However, 
in the case of balloon-borne detectors, the uncertainties 
introduced by large corrections for secondary electrons 
produced m the residual atmosphere above the instrument 
may contribute significantly. The satellite experiments 
(21) (45) have different, but no less severe, problems 

correcting for background events. In our work we make use 
of new calculations of the growth curves of atmospheric 
secondary positrons and negatrons (47) , together with data 
gathered during ascent, to determine the atmospheric 
secondary contribution at float altitudes. (A complete 
discussion is found in Chapter IV, Section E) . We derive 
an extraterrestrial electron spectrum which is in general 
close to or somewhat below previously reported intensities. 
Our results connect well to the most reliable data measured 
below and above | this energy interval. Our measured extra- 
terrestrial positron and negatron spectra give a charge 
ratio which is substantially constant throughout the 
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mterval and indicate that the collision source contributes 

i 

30-50 percent of the measured electron flux over most of 
the range between 11 and 204 MeV. We derive the absolute 
solar modulation of positrons by comparing our measured 
spectrum with the interstellar positron spectrum calculated 
for the collision source. Neglecting energy loss during 
modulation we derive the interstellar electron spectrum 
between 11 and 204 MeV and extrapolate this spectrum to 
higher energies. This result is compared with the electron 
spectrum inferred by other investigators from radio data. 
The electron spectrum from primary sources alone is also 
derived. We examine the effect of adiabatic deceleration 
of cosmic rays passing through the solar wind and find 
that many of our conclusions regarding the interstellar 
electron spectrum remain essentially unaffected .even for 
relatively large energy loss. 



II. INSTRUMENT 


A) Detector system 
1) Overall description 


The detector is a magnetic spectrometer, an instru- 
ment which measures the charge sign and magnetic rigidity 
(momentum divided by charge) of charged particles by 
determining their deflection in a magnetic field. Instru- 
ment parameters are chosen such that the rigidity interval 
6-200 MV can be studied with good resolution. A schematic 
cross-section of the detector is shown in Figure 1. 

A triple coincidence of Telescope Counter #1 (Tl) , 

v 

Telescope Counter #2 (T2) , and the Lucite Cerenkov Counter 
(C) , triggers the spark chamber high voltage and initiates 

i 

the data readout cycle. Two 4-gap wire spark chambers de- 
termine the trajectory of the particle before and after 
traversal of the gap of a permanent magnet. The spark 
locations are detected by magnetostrictive pickups and re- 
corded in digital form on magnetic tape. 

v 

The Lucite Cerenkov counter is designed to discrimi- 
nate against upward-moving splash albedo particles and to 
eliminate the nucleonic component of the cosmic rays in the 
rigidity interval of interest. The velocity threshold for 

V j 

Cerenkov emission in lucite (index of refraction = 1.49) 
is 0.67 c; corresponding rigidity thresholds are .46, 845 
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TOP GUARD 



Fig. 1 . Cross-section of the detector system. 
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and 1690 MV for electrons, protons, and alpha partrcles, 
respectively. Electronic pulse discrimination increases 
the rigidity threshold an additional 15 percent. 

The specially designed Magnet Guard Counter, MA, is 
shown in Figure 2, together with the magnet, in an "explo- 
ded" isometric projection. The magnet guard counter covers 
the top of the magnet and completely lines the gap volume 
leaving a 3-cm x 12-cm open passage. This guard counter 
eliminates particles which might interact or scatter in the 
magnet pole pieces and, together with T1 and T2, defines 
the acceptance cone of the detector. Additional guard 
counters surround the sides of the instrument and cover the 
top with the exception of the telescope aperture. These 
counters eliminate charged particles which enter the detec- 
tor from outside the acceptance cone and might subsequently 
interact, producing particles which trigger the telescope 
counters. All guard counters are m active anti-coincidence. 

The atmospheric pressure during flight is monitored 
by a Wallace-Tiernan aneroid barometer, (FA 160) which is 
photographed, together with a clock and a thermometer, at 
5-minute intervals. The calibrated barometer is accurate 
to within ± 0.1 mb at 2.4 mb, our typical float altitude. 

- An additional low-pressure gauge (Metro-physics Inc, Santa 
Barbara, Calif) is read out electronically. 
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Fig. 2. Exploded view of the magnet. 


and magnet guard counter. 
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The entire instrument, with the exception of the 
barometer and 30 silver-zinc batteries, is enclosed in a 
pressure-tight aluminum gondola during flight, A 0,012- 
inch| aluminum window allows particles incident from above 
to penetrate the gondola with minimum scattering. The gon- 
dola is insulated and the temperature within is maintained 
above 45°F by thermostatically-controlled heaters. T£e 
weight of the instrument package is 370 lbs. 

2) Magnet 

The analyzing magnet (Indiana General Corp., Valpar- 
aiso, Ind.), which is shown in Figure 2, has a 4-cm air gap 
between 9-cm x l4-cm rectangular pole faces. Alnico 8 
permanent magnets are used in a magnet circuit designed to 
minimize external stray fields. At the position of the 
photomultiplier tube for the magnet guard counter, for 
example, the field is less than one gauss. The nominal 
flux density within the magnet gap is 1000 gauss. The 
three orthogonal components of the field were measured at 
one-cm intervals throughout the volume accessible to 
particles out to a distance 5 cm above and below the magnet 
Beyond these joints readings were taken at larger intervals 
The accuracy of individual field measurements is estimated 
at about 3 percent, based on the scatter of adjacent 
readings from a smooth curve and on deviations from 
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expected symmetries. Figure 3 shows the flux density 
measured along several vertical paths. Locations are 
specified in a righthanded coordinate system with the ori- 
gin at the center of the magnet gap. The z-axis is vertical 
and positive upward; the x-axis is perpendicular to the 
pole faces and positive toward the south pole. (This coor- 
dinate system is used throughout this thesis'. It is shown 
explicitly in Figure 5) . A permanently mounted Hall effect 
device (F. W. Bell, Inc, Columbus, Ohio) was monitored as 
part of the check-out procedure before and after each 
flight; no change m the field strength was ever noted. 

The geometrical factor, G, of the detector varies 
with rigidity due to the effect of the magnetic field on 
the beam. A Monte Carlo-type calculation determined G at 
six different rigidities R from 6 to 200 MV. A trajectory 
was selected at random from a simulated uniform isotropic 
flux of particles of rigidity R^ incident on counter Tl. 

The trajectory was traced through the detector utilizing 
the measured magnetic field; it was rejected if it failed 
to pass through counter T2 or the aperture in the magnet 
guard counter. This procedure was repeated until 1000 
trajectories had been accepted. The geometrical factor 
at R^ is then 
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Fig. 3. Magnet flux density an the gap of the analyzing 
magnet vs. position. The curves represent the 
field components along three vertical paths through 
the magnet gap. See Figure 5 and the text on page 
14 for a description of the coordinate system 
employed. 


1. Solid curve: 

x = 

0 

y - 0.5 

cm 

2. Dashed curve: 

x = 

0 

y = 5.5 

cm 

3. Dot-dashed curve: 

x = 

1 cm 

y = 5.5 

cm 


B z and By for paths 1 and 2 are identically zero 
and are therefore not shown explicitly in the 
figure . 



FLUX DENSITY (GAUSS) 


l 


- 16 - 



Z-COORDI NATE (CM) 
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1000 

G(R l ) = N(R i ) A “ (1) 

where A is the area of Tl and S3 is the solid angle out of 
which the N (R^) initial trajectories were chosen. The 
geometrical factor as a function of rigidity is shown m 
Figure 4. The calculation is accurate to about 3 percent 
and within this accuracy the geometrical factor is the same 
for both positive and negative particles. 

For the following discussion reference should be 
made to Figure 5 where the relevant trajectory parameters 
are shown. We define entrance and exit parameters of the 
trajectory at convenient planes z — ±z Q/ above and below 
the region of field. The angle ri is measured between the 
trajectory and the y-z plane; ai and a 2 are the entrance 
and exit angles projected in the y-z plane and measured 

f 

from the negative z-axis. The deflection angle 0 is defined 

m the y-z plane and is given by 
I 

0 = a - a, (2) 

2 i 

Angles measured clockwise are taken to be positive. Since 
the magnetic field is in the positive x-direction, the 
deflection angle is then positive for positively charged 
particles traversing the field in the downward direction. 



GEOMETRICAL FACTOR (CM -SR) 



Geometrical factor of the detector vs. particle rigidity. 


Fig. 4. 




' Fig. 5. Schematic view of a particle trajectory through the magnet gap, showing 
reference coordinate system and trajectory parameters used in the 
calculation of the particle rigidity. 
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The deflection of a particle of rigidity R moving 
through a magnetic field B is given by 

0 = 3xl0 ~ 4 / B. dl (3) 

R ± 

where Bj_ is the component of B normal to the trajectory and 

dl is an increment of distance. The deflection is given in 

radians for R, B, and 1 in MV, gauss, and cm, respectively. 

The line integral is evaluated along the trajectory of the 

particle and is known as the magnetic path, M. If the 

field and the trajectory of the particle are known, M can 

be evaluated and eq. (3) solved for the rigidity. In 

general M will differ along different possible trajectories 

through the magnetic field. In Figure 6 we show the range 

of values of R times 0 for the random trajectories accepted 

m the calculation of the geometrical factor. Plotted are 

the mean value, the rms deviation (solid bar) , and the 

! 

extreme values (dashed bar) for each of the sample popula- 

j 

tions in that calculation (1000 trajectories for each 
rigidity) . At all points the mean of R x 0 is within one 
percent of 3.55 MV; the rms deviation is typically 2 per- 
cent of the mean and contains ~ 70 percent of the sample 
while the extreme values differ from the mean by about 7 
percent. The momentum resolution of the detector is 
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Fig. 6. Range of trajectory deflection angles in the 
detector as a function of particle rigidity. 
The mean value, rms deviation (solid bar), 
and extreme values (dashed bar) of rigidity 
x deflection angle are shown for a random 
distribution of 1000 incident trajectories 
at each of 6 rigidities. 



ANS) 
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^ 30 percent, FWHM (See Chapter IV, section A) . Within 
this resolution we can approximate eg. (3) by 

R = mv (4) 

0 


with negligible error. 

Equation (4) has no dependence on the angle n . 
Furthermore, the complete sheathing of the magnet by the 
magnet guard counter eliminates the necessity of protecting 
against ‘interactions in the pole faces by detailed trajec- 
tory reconstruc'tion. This permit's us to read out spark 
locations in the y-z projection only, saving both data 
storage and detector live time. 

3) Scintillation Counters 

i 

The two telescope counters Tl and T2 and the nine 
guard counters are constructed from NE 102 plastic scin- 

i 

1 

tillator (Nuclear Equipment Corp., San Carlos, Calif.) and 
employ RCA 4439 photomultiplier tubes (RCA, Harrison, N. J.). 
Except for the specially cast magnet guard counter, all of 
the counters are flat sheets. Individual dimensions are 


listed in Table 1. 



TABLE 1 


SCINTILLATION COUNTERS 


Desig- 






- 

Connection to , 

nation 

Function 


Dimensions 

(Inches) 

photomultiplier 

__T1 

Telescope 1 

1/8 

X 

4 

x 7 


Lucite light pipe 

T2 

Telescope 2 

1/8 

X 

4 

x 6 (a) 


Lucite light pipe 

Top 1 

Top Guard 

3/8 

X 

15 

1/2 x 

17 1/8 (b) 

Direct coupling of 
photomultiplier to 

Top 2 


with 

hole 5 

x ,8 3/8 

scintillator for all 








guard counters 

Side 1 


3/8 

X 

10 

1/2 x 

15 1/4 

1 

to 

Side 2 


3/8 

X 

10 

1/2 x 

15 1/4 

4^ 

1 

Side 3 

Side Guard 

3/8 

X 

10 

1/2 x 

15 1/4 


Side 4 


3/8 

X 

10 

1/2 x 

15 1/4 


Side 5 


3/8 

X 

15 

7/16 x 22 5/8 


Side 6 


3/8 

X 

15 

-7/16 x 22 5/8 


MA 

'Magnet Guard 

‘ (c) 







(a) T2 is not rectangular: it is a 4-inch wide strip cut from the center of a 6-inch 

diameter disc. 


(b) Dimensions are for entire top guard; there are two "L" 

(c) Specially cast; see Figure 2 and text on p. 25. 


-shaped pieces. 
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The photomultiplier tubes for counters T1 and T2 
are mounted to lucite light pipes which view the scintil- 
lators edge-on. The resolution for minimum ionizing 
particles is 60 percent, FWHM, for both counters. The 
photomultiplier of each guard counter xs coupled directly 
to the large flat side of the scintillator near one corner. 
Optical coupling is made with Dow-Corning 20-057 optical 
grease (Dow-Corning Corp., Midlands, Mich.). Despite the 
large size of the exterior guard counters the light collec- 
tion efficiency is relatively independent of the point of 
passage of the particle outside a circle of •~ r 10-cm radius 
centered on the phototube, and is everywhere sufficient to 
allow conservative discrimination levels to be set. 

The magnet guard counter, which is shown m Figure 2, 
is cast in one piece m order to minimize light loss. The 
scintillator is 9.4-mm thick except for the area directly 
covering the pole faces which is 5-mm thick. The upper 
edge of the gap lining is beveled to aid reflection around 
the corner. Resolution is 90 percent for perpendicular 

s 

traversal of a single thickness of the pole face sheath by 
a minimum ionizing particle. The discrimination level is 
set to trigger on 99.8 percent of such particles. Particles 
which trigger counters Tl, T2, and C and have interacted or 
scattered m the magnet material must have passed either 
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through the top surface or obliquely through the magnet 
gap lining of MA. In either case, the light output is sig- 
nificantly greater than for perpendicular traversal of the 
pole face sheath. We therefore assume total rejection of 
such events. 

The presence of the magnet guard counter as one of 
the elements defining the acceptance cone produces an 
uncertain rejection zone at the edges. This is due to the 
finite path length in MA which is required for rejection of 
the event. This zone is about 0. 5-1.0 percent of the total 
geometrical factor of the instrument. While negligible as 
an error in the geometrical factor, it represents an addi- 
tional source of scattering for 0. 5-1.0 percent of the 
'accepted electrons. The effect is included in the angular 
scattering distribution discussed in Chapter IV, section A. 

4) Cerenkov Counter 

1 

The Cerenkov radiator is a 1-inch thick, 7 5/16-inch 
diameter disc of ultraviolet-transmitting lucite. It is 
mounted directly on the face of an EMI 9623 B photomulti- 

s 

pier tube (Whittaker Corp., Plainview, N.Y.). This tube is 
7 1/2 inches m diameter and has a quartz window. Optical 
coupling is made with Sylgard 182 resin (Dow-Corning, 
Midland, Mich. ) . 
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Resolution of the Cerenkov counter is 45 percent, 

FWHM. The discrimination level is set at 25 percent of the 
output from relativistic muons. The detection efficiency 
for relativistic particles moving backward through the 
Cerenkov counter is ~4 percent. During pre- and post- 
flight checkout, the singles rate of the Cerenkov counter 
was always <300 counts/sec and probably did not exceed 
~1000 counts/sec during flight. Since the resolving time of 
the triple coincidence is 1 ysec, during flight the pro- 
bability was less than ~10 that a random pulse in C' 
would accompany a T1 T2 coincidence. 

Even though electrons from 6 to 200 MeV are well 

V 

above the Cerenkov threshold, at the lowest energies the 
output and resolution of the counter are degraded somewhat 
due to increased scattering of the electrons. There is a 
higher probability of backscattering or stopping in the 
radiator which result in a shorter effective path for 

V 

Cerenkov radiation directed toward the phototube cathode. 
There are no directly applicable experimental data which 

j 

allow us to determine the effect on our detection efficiency. 
In the work of Jupiter, Loner gan, and Merkel (48) with 8.0 
MeV electrons incident on a '2.4 g/cm aluminum absorber 
(0.1 radiation length), the backscattered flux and the 
flux leaving th'e absorber in the forward hemisphere are 
respectively 2.7 percent and 77.3 percent of the incident 
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flux. The other 20 percent have a path length which, be- 
cause of scattering, is sufficiently long that the particles 
lose all of their energy by ionization loss. Our Cerenkov 
radiator is one inch of lucite, which is 2.6 g/cm^ or 0.06 
radiation lengths. With our threshold set at 25 percent of 
the muon output, a radiating particle need only pass through 
^1/4" of lucite in order to have 50 percent probability of 
detection. Scaling from the data of Lonergan et al . we 
estimate that all but about 5 percent of the 8 MeV electrons 
will penetrate 1/4" (.65 g/cm 2 ) of the lucite. In addition, 

a scattering electron will have, in general, a longer effec- 

v 

tive path for radiation of Cerenkov light than the depth 
to which it actually penetrates. It follows that the de- 
tection efficiency for 8.0 MeV electrons is probably 95 
percent. Because the exact value is uncertain but appears 
to be greater than 95 percent even at this unfavorable 
energy, we have assumed the detection efficiency to be 100 
percent over the entire interval from 6 to 200 MeV. 

5 ) Spark Chambers 

The two identical wire spark chambers each consist 
of 4 independent and self-contained modular gaps. Figure 
7 shows an exploded view of a single module. The three 
frame pieces are machined from Supramica 500, a type of 
bonded mica (Mycalex Corporation of America, Clifton, N.J.). 
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Fig. 7. Exploded view of a spark chamber module. 
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The upper and lower frames are each 1/8" thick and serve 
as insulating, spacers between adjacent gaps. The central 
frame determines the spacing of the wire planes, which is 
.250", and the active area, which is 5" x 9". The high 
voltage and ground planes employ .0022" diameter silver- 
coated beryllium copper wire (Little Falls Alloys, Pat- 
terson, N.J.). The wires are parallel and evenly spaced 
48 per inch and are held in place with Shell 828 epoxy 
(Shell Oil Co., New York, N.Y.}. The 5 edgemost wires 
within the active area at either end are raised slightly 
and are looped by a strip of aluminized mylar which extends 
over the inner edge of the central frame. This feature 
prevents a high field concentration near the edge which 
could lead to spurious sparking. 

All the wires of a single plane are connected at one 
end to a buss bar to which external connection can be made. 
Under the ground plane wires between the buss bar and the 
active chamber area is a slot which accepts the magneto- 
strictive pickup. (Use of a magnetostrictive technique for 
spark chamber readout has been described by Perez-Mendez 
and Pfab (49)-). As the pickup we use a .004" x .009" 
ribbon of the magnetostrictive material Reraendur (Wilbur B. 
Driver Co., Newark, N.J.) which is held near the wires of 
the ground plane by an aluminum "wand". The wand slips 
into the module frame as an independent sub-assembly. The 
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wand and ribbon are insulated from the chamber wires by a 
.001" layer of mylar. The current pulse from a spark m 
the chamber passes over the wand inducing a mechanical 
deformation in the ribbon. The mechanical pulse travels 
away from the point of origin at the speed of sound in 
Remendur, which is 5.3 cm/ysec at room temperature. At 
each end of the ribbon rubber pads damp the pulses, pre- 
venting reflections. Near one end the ribbon threads a 
small 200-turn pickup coil. Small bias magnets produce a 
magnetic field m the ribbon at the coil. Passage of the 
mechanical pulse induces a voltage signal in the coil by 
the inverse magnetos trictive effect. 

Two fiducial wires, which carry pulses each time the 
high voltage is triggered, pass over the wand outside the 
chamber wires, one near either end. The quantity which is 
measured is the time delay between the arrival at the pick- 
up coil of the nearer fiducial pulse and the subsequent 
spark pulse. Since the velocity of propagation of the 
pulses is known, the distance between the first fiducial 
wire and the wire carrying the spark current may be easily 
calculated. With this technique all spark location measure- 
ments are made relative to a standard fixed in the chamber. 
Furthermore, by separating the pickup coil and the nearer 
fiducial wire sufficiently, the fiducial pulse will not 
arrive at the coil until all RF noise from the spark has 
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disappeared, and hence the timing circuitry need not be 
noise -immune . The presence of the second fiducial provides 
a known constant distance with which to monitor changes of 
the propagation velocity due to temperature or other effects. 
It also performs a check of the timing and digitizing 
circuitry. 

After all wires have been attached and all electrical 
connections made, the upper and lower frame pieces are 
epoxied to the central piece. The four modules of each 
chamber are mounted in pairs, with their high voltage 
planes back to back, on the top and bottom of a 4.25"-high 
open aluminum box. The gas volume is closed off by 0.5-mil 
mylar - 0.5-mil aluminum laminated foils on the top and 
bottom of the chamber. The gas used is standard "spark 
chamber neon", which is 90 percent neon and 10 percent 
helium. A 2 percent ethanol admixture acts as a quenching 
agent. During flight gas flows continuously through the 
chambers at a rate of 0.25 ml/sec and is exhausted into the 
gondola volume. 

The typical efficiency for detection of a single 
charged particle is 96 to 98 percent in each gap. During 
some flights a few of the modules developed edge sparks 
which did not, however, seem to affect the efficiency of 
sparking at the! trajectory location. The distribution of 
the deviation of the measured spark location from the true 
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trajectory is approximately gaussian with a standard 
deviation 


- .008" (5) 

This excellent spatial resolution, which is less than the 
wire spacing, is the result of the magnetostructive read- 
out, which automatically indicates a mean location, weighted 
by the relative currents, should two or more wires carry 
the spark current. 

The data readout determines the wires which carry 
spark current. It follows that a consistently undistorted 
projection of the trajectory requires that the wires of all 
modules be parallel. The accuracy of the wire alignment is 
on the order of ± .3 milliradians within a single chamber 
and about 1 milliradian between the two spark chamber 
boxes. The spark chambers and the magnet are mounted 
rigidly together so that no shifting can occur after the 
alignment is made. 

B) Electronic system 
1) Overall description 

Figure 8 shows the general block diagram of the 
electronics system. Data collection is broken up into 16- 
mmute cycles each of which consists of a 15-minute segment 
called Phase A followed by a 1-minute segment called Phase B. 
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During phase A, a triple coincidence of counters Tl, 

/ 

T2, and C, unaccompanied by a pulse in any of the guard 
counters, triggers the high voltage pulsers and initiates 
the data readout cycle. A busy signal blocks the coinci- 
dence during the 350 milliseconds required to process the 
event. During the first 12 ysec all noise sensitive cir- 
cuitry is held in the reset state. During the subsequent 
1000 usee the spark locations, atmospheric pressure, and 
gondola temperature data are digitized. Only the location' 
of the spark nearest the first fiducial is stored. The 
presence of subsequent sparks (not including the far fidu- 
cial pulse) is indicated by a "Multiple Spark Indicator" 

(MSI) bit. These data, together with the time information, 
are then stored on 16-channel magnetic tape. The tape 
recorder runs continuously at 10 inches per minute and holds 
sufficient tape for a 30-hour flight. During ground 

* 1 I 

testing the data can also be recorded directly on punched 
paper tape, which allows manual readout. 

During phase B the normal coincidence trigger input 
is blocked. After sufficient time is allowed to finish 
recording any preceding event, a trigger pulse is initiated 
internally. Since this trigger is not associated with any 
real particle in the detector, no sparks should occur in 
the chambers arid the digitized location should be that of 
the far fiducial. In practice a spark chamber gap sometimes 
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breaks down in fandom locations, but the fiducial appears 
frequently enough to provide a check of the digitizing 
operation. 

The time remaining in phase B is used to monitor 
counter rates. There are 4 rate scalers: Tl T2 double 

coincidence rate, Tl ~ T2 A C triple coincidence rate, 
magnet guard counter singles rate, and total guard counter 
singles rate. Throughout phase A these scalers are 
blocked. At the beginning of phase B they are reset and 
then count until the end of the one minute period, at which 
time the accumulated count is recorded on the magnetic 
tape. The rates are monitored in order to check possible 
counter malfunction, to detect variations in background 
radiation which might be caused, for example, by solar 
flares, and to determine detector dead time due to the 
guard counters. By blocking the spark chamber trigger input, 
the necessity for shielding the scalers from spark generated 
RF noise is eliminated. Only the clock and time scaler 
must be shielded, which is simple since no fast signals are 
involved. The clock, which provides the timing for the 
Phase A-Phase B cycle, is a Bulova Accutron. 

Total power consumption of the electronics system is 
fourteen watts. When the heaters are all on an additional 
80 watts is dissipated in the gondola. Power is supplied 
by 30 silver-zinc batteries (Yardney Electric Corp., Los 
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Angeles, Calif.) and 4 nickel-cadmium batteries (Gould, St. 
Paul, Minn.), The battery packs provide voltages between 
-6v and +30v. High voltages for the photomultiplier tubes 
and spark chambers are produced by DC-DC converters 
(Crestromcs , Crestline, Calif., and Mil Associates, Hudson, 
N.H.) . 

2) High voltage pulsers 

A 4.5 kilovolt pulse is applied to the spark chamber 
by pulsers utilizing KN 22 krytron tubes (E. G. & G. , Inc., 
Salem, Mass.). There are eight pulsers, each one of which 
drives a single module and its associated fiducial wires. 

The pulser circuit and network are shown in Figure 9. The 
resistor R in parallel with the spark chamber module CM is 
used to adjust the decay time of the applied pulse for 
optimum performance. Too low a value causes sparking 
efficiency to decline while too large a value produces in- 
creased spurious breakdown. Typical values vary between 
200 and 600 fi. The value of the capacitor resistor chain 
in the fiducial lines is chosen to give magnetostrictive 
pulses .approximately equal in amplitude to average spark 
discharge pulses. 

The pulsers are all triggered by a single avalanche 
transistor-puls, e transformer circuit which is itself 
directly triggered by the output of the coincidence circuit. 
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Fig. 9. High voltage pulse network. 

CM = spark chamber module. 
R = adjustable resistor. 
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fiducial wires. 
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The 10-90 percent risetime of the HV pulse at the spark 
chamber is 30 nsec. The total delay from passage of the 
particle to appearance of the high voltage pulse is- about 
275 nsec which is divided as follows: 


Photomultiplier Tubes 
Coincidence Circuitry 
Avalanche Circuit 
Pulser 


100 nsec 
60 nsec 
80 nsec 

35 nsec 
275 nsec 


3) Spark chamber readout 

The output of the magnetostrictive readout pickup 
coil is a sequence of pulses corresponding to the two fidu- 
cials and whatever sparks may have occurred in the spark 
chamber gap. These pulses have the shape shown in Figure 
10A(a). Despite being relatively broad (the central lobe 
is about 400 nsec wide, equivalent to 2 mm on the magneto- 
strictive ribbon) . The symmetry and reproducibility of the 
pulse shape allow the separation between pulses to be- 
determined to better than 50 nsec. Our method is similar 
to that of Kirsten, Lee, and Conragan (50) . After the 
small signals from the coil are amplified by a preamplifier 
mounted directly on the wand, they go to the so-called 
"magnetostrictive discriminator." Here the sequence of 
pulse transformations shown in Figure 10A takes place. 




(b) 



(c) 



(d) 
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B. Block diagram of spark chamber readout electronics. 


Fig. 10. Magnetostrictive spark chamber readout 
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The input pulse, (a) , first has its outer lobes clipped, 

(b) , after which it is differentiated, (c) . The output 
(d) of the zero-crossing discriminator which follows is a 

i 

pulse with a relatively slow leading edge (100-150 nsec 
rise time) but a fast trailing edge (10 nsec fall time) . 

The trailing edge corresponds to the zero crossing of (c) 
and hence to the peak of the input pulse (a) . This edge 
is used in the subsequent routing and digitizing elec- 
tronics, a block diagram of which is shown in Figure 10B. 
This part of the circuitry is similar to the commercially 
available Model 180 Multiple Time Digitizer (LeCroy 
Research Systems Corp., Elmsford, N.Y.) but has been 
adapted to our special requirements. Each module feeds 
into a single 11-bit scaler. Initially the scaler is 
reset and gates 1 and 2 are blocked. The arrival of the 
first pulse (the near fiducial) opens gate 1 which allows 
pulses from the 20 MHz oscillator to reach the scaler. 
Scaling continues until the arrival of a second pulse 
again blocks gate 1. A third pulse opens gate 2 and a 
fourth sets the multiple spark indicator. Any subsequent 
pulses have no effect. The indicator bit is thus set 
only if 2 or more pulses occur in the gap. The external 
enable allows digitization only in the 75 jusec immediately 
following an event. This protects against random noise 
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pulses setting the multiple spark indicator during the 
long (240 msec) period required to record the event. 

The spark information recorded for each module is 

1. the time delay between the first and second 
pulses on the magnetostrictive ribbon, digitized 
to 50 nsec (equivalent to .25 mm), and 

2. an indication if more than one spark occurred 
in the gap. 



III. BALLOON FLIGHTS 


The data reported in this thesis were derived from 
three high-altitude balloon flights launched from Fort 
Churchill, Manitoba, on July 16, July 21, and July 29, 1968 
(Universal Time) . The flights are designated Cl, C2, and 
C3, respectively. Relevant flight information is summarized 
in Table 2. 

Figure 11 shows the trajectories of the three 
flights. Invariant latitude contours, calculated from 
the internal field only (51) , are included to show the 
trajectories in the geomagnetic field. 

In Figure 12 we show the altitude curves for the 

flights. In each case the launch was timed so that the 

2 

detector would pass through the level at 100 g/cm residual 
atmosphere after the evening transition to the low night- 
time geomagnetic cutoff (see Chapter IV, section E and 
Chapter V, section A) ; the evening transition occurs at 
about 18:00 local time (00:20 UT) . The float altitude 
during the nighttime interval stayed within the range 
2. 2-2.6 g/cm 2 residual atmosphere for all three flights. 

In Figure 13 we show the daily average of the 
hourly count rate for both the Churchill (52) and the 
Mount Washington (53) neutron monitors, as well as the 
daily average of the 3-hour range indices for the 



TABLE 2 


BALLOON FLIGHTS 


Flight Number 

Cl 

C2 

C3 

Launch Date (1967) (a) 

July 16 

July 21 

July 29 

Launch Time (a) 

00:27 

01:04 

01:31 

Reach 100 g/cm 2 (a) 

01: 32 

02:00 

02:42 

Begin Float (a) 

03:49 

04:09 

04:53 

Terminate (a) 

17:15 

19:18 

01:11 (July 

Floating Depth (g/cm 2 ) (b) 

2. 2-2. 6 

2. 2-2. 6 

2. 1-3. 4 

Total Sensitive Time at Float (Min) 

625 

695 

932 

End Night Interval (a) (c) 

09:09 

09:45 

10:45 

Average .Floating 

Depth — Night (g/cm 2 ) (c) 

2.45 

2.5 

2.35 

Total Sensitive Time at 

Float — Night (Min) (c) 

251 

254 

265 

K (d) Total Flight 

p Nighttime Only 

2+ 

' 2- 

1+ 

2 

1 

' 1+ 



Table 2 Continued 


Mt. Washington Neutron 


Monitor (e) 

Total Flight 

2163 

2200 

2227 


Nighttime Only 

2147 

2194 

2230 

Churchill Neutron 

Monitor (e) 

Total Flight 

5857 

5942 

6025 


Nighttime Only 

5861 

5921 

6003 


(a) Universal Time 

(b) See Figure 12 

(c) For the period included in the nighttime data. For actual transition times see «£• 

cn 

Chapter V, section A. 

(d) Mean of three-hour indices during float (54) 

(e) Mean of hourly count rate during float (52) (53) 




FLIGHT Cl 

FLIGHT C2 

FLIGHT C3 


Fig. 11, Trajectories of the balloon flights. Invariant 
latitude contours are derived from the internal 
geomagnetic field only (Ref. 51). 
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Fig. 12. Altitude curves. The curves are labeled 
with the flight number. 



m 


1 


u x 

O LU 

< Q 

or 


< x 


9* I 


6100 


6000 


■j 


5900 


5800 


CHURCHILL NEUTRON 
MONITOR 

C3 

C2 ^ 


2200 


2150 


MT. WASHINGTON 
NEUTRON MONITOR 


21001 


JULY, 1968 

Pig. 13. Planetary magnetic index, Kp, and Churchill 

and Mt. Washington neutron monitor count rates 
during July, 1968. The numbered arrows 
indicate the times of our flights. 
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planetary magnetic index, K (54), for July 1968. All of 

P 

our flights occurred during or just after the recovery 
phase following a Forbush decrease. The average hourly 
count rate of the Churchill neutron monitor was 5857, 5942, 
and 6025 for flights Cl, C2, and C3, respectively. These 
rates differ among themselves by about 3.2 percent. For 
comparison, the average neutron monitor count rate during 
July 1967 was about 4 percent higher than the average rate 
during July 1968. 

The solar proton monitor aboard the Explorer 34 

satellite recorded no flux above 30 MeV and only negli- 

o 

gible flux (0.2 particles/cm *360*51) above 10 MeV 
throughout flights Cl and C3 (55) . The same is true for 
the earlier part of C2, but at 15:00 UT a large short- 
term increase is evident. The event produced a significant 
f luxv of protons above 60 MeV. Activity is also evident 

in the K indices at about the same time and a principal 
P 

magnetic storm occurred at about 16:00 UT (52). The 
relationship to our measured data is discussed in Chapter 
V, section A. 



IV. DATA ANALYSIS 


A) Detector Resolution 

Several effects contribute in determining the mag- 
netic rigidity resolution of our detector. The results 
of limiting the trajectory measurement to a single plane 
projection and of using a constant magnetic path approxima- 
tion have already been discussed in Chapter II. An uncer- 
tainty of a few percent is introduced m the rigidity 
determination by these simplifications. More significant 
uncertainties at all rigidities are the result of the in- 
trinsic angular resolution of the instrument and/or multi- 
ple electron scattering. The intrinsic angular resolution 
determines the ability of the detector to measure the 
actual deflection of the particle trajectory. Multiple 
electron scattering adds a random angular deviation to the 
deflection due to the magnetic field. It therefore deter- 
mines how greatly the trajectory deflection angle differs 

I 

from the angle corresponding to the 'particle rigidity R 
according to eg. (4), Chapter II. We shall discuss the 
errors in terms of the deflection angle 0, since the error 
distributions are symmetric in angle. Because the 
rigidity R is inversely proportional to 0, the equivalent 
distributions in R are skewed somewhat toward higher 
rigidities. Each of the two effects will be discussed 
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separately. Both will then be combined to give the "true" 
angular resolution, which corresponds to the ability of 
the detector to' measure the rigidity of a particle. 

The discussion in this section and the next uses a 
number of distributions for which standard deviations are 
defined. Since it has not been possible to employ self- 
evident symbols in every case, we have summarized in Table 
3 the symbols used for the various standard deviations, 
together with a short definition for each and the number 
of the equation where the symbol is introduced. More 
complete definitions are found in the text. 

We consider first the error in the measured deflec- 
tion angle due to the intrinsic angular resolution of the 
detector, assuming no scattering. The intrinsic angular 
resolution defends on the ability of the instrument to 

i 

define the entrance and exit angles, a-^ and a 2 / °f the 
particle trajectory (see Figure 5). It is independent of 
the rigidity or the species of the particle. Within each 
spark chamber the trajectory is determined by making a 
least-squares fit of the measured spark locations to a 
straight line. With our spark plane arrangement (see 

< \ i 

V • 

Figure ,1), the standard deviation a c of the determination 
of the slope (tan a) of the particle trajectory m a single 
spark chamber is given approximately by 
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TABLE 3 

SUMMARY OF SYMBOLS USED FOR STANDARD DEVIATIONS 

Symbol Definition Where first 

used (a) 

<?£ Standard deviation (S.D.) Eg. (5), Ch. II 

of the measured spark 
location in a single spark 
plane 

u a S.D. of the measured entrance Eq. (3) 

and exit angles, and 

u g S.D. of the measured entrance Eq. (1) 

and exit slopes, tan and 
tan (*2 

cr A S.D. of the measured deflec- Eg. (4) 

tion angle arising solely 
from the intrinsic angular 
resolution of the detector 

S.D. of the calculated pro- Eg. (5) 

jected scattering-angle 
distribution for ielectrons 

Uq S.D. of the measured deflec- Eq. (7) 

tion angle 0, including 
electron scattering and the 
intrinsic angular resolution 
of the detector 



S.D. of the 
(see Figure 

measured 

18) 

angle w 

Eq. 

(16) 

<*A 

S.D. of the 
(defined in 

measured 
Eg. (13) 

angle A 

Eq. 

(16) 


(a) Equation numbers refer to Chapter IV unless otherwise 
indicated. 

| 
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for good sparks in 4 planes 

( 1 ) 

for good sparks m 3 planes 

where 0 ^ is the standard deviation of the distribution of 
the measured spark location about the true trajectory 
position and r is the overall separation within the spark 
chamber. In this case r = 5.3" and 0 - .008" (eq. 5, 

A> 

Chapter II) and hence 



{ . 0015 for 4 planes 

:0018 for 3 planes 


( 2 ) 


The relationship of 0 to the standard deviation 0 ^, for the 
incident angle a depends weakly on the value of a within 
the acceptance cone of the instrument. We shall use a 
value for a a obtained by averaging over the possible tra- 
jectories within the acceptance cone. The result/ valid 
for both the entrance and exit angles, ctj and is 


0 


a 


.0019 radians for 4 planes 

- 

.0022 radians for 3 planes 


(3) 
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The error distribution for a will be approximately gaussian, 
since the distribution of the measured spark location for 
any spark plane is approximately gaussian. - The deflection 
angle 0 is the difference of and (eg. (2) , Chapter 
II) . It follows that the standard deviation u A for the 
error m the measurement of the deflection angle due solely 
to the intrinsic angular resolution of the detector is 

.0027 radians for 8 planes 
.0029 radians for 7 planes (4) 

.0031 radians for 6 (3+3) planes. 

For data analysis, events are divided into 3 
categories : 

' 1. Those events where all 8 spark planes define the 
trajectory, referred to as "perfect" events; 

i 

2. Those events where a single spark plane m one 
or both spark chambers malfunctioned (registering either no 
spark or a spark off the trajectory defined by the other 3 

I 

planes), referred to as "one-error" events; and 

3. Those events where more than one spark plane in 
one or both spark chambers malfunctioned, referred to as 
"multiple-error" events. Only the first two categories are 
used for the determination of particle rigidity. (See 
section B for more details on event selection criteria.) 


“a ' / 


ct a 1 +a a 2 


a < 
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In Figure 14 we show the measured distribution of 
deflection angles for a mono-energetic beam of positrons of 
800 MeV from the California Institute of Technology 
Synchrotron. At this energy the deflection resolution is 
determined primarily by the intrinsic angular resolution of 
the detector. Also shown is a gaussian distribution for 
which 0 = .0028 radians. The agreement of the measured 
distribution with eg. (4) is seen to be quite good for both 
"perfect" and "one-error" events. 

We shall consider now the effect of multiple scat- 
tering, ignoring the intrinsic angular resolution of the 
detector. The scattering which affects the deflection 
measurement occurs principally at the wire spark chamber 
planes and foil covers immediately above and below the 
magnet. In Figure 15 we show (curve 1) the distribution of 
the projected scattering angle <f> calculated for electrons 
of momentum p according to the theory of Moliere (data 
from ref. 56). Both the differential distribution, dN/d4>, 
normalized by dividing by p, and the integral distribution, 
N (><{>}, are plotted versus pcj>. Plotted in this way curve 1 
can be used for all electron momenta above a few MeV/c. It 
follows that if k is the standard deviation of the distri- 

f 

bution as plotted, then the standard deviation of the 
projected scattering-angle distribution is 
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DEFLECTION (MILLIRADIANS) 


Fig. 14. Distribution of measured deflections in an 
800 MeV positron beam. "Perfect" events 
(solid histogram) and "one-error" events 
(dashed histogram) are plotted separately. 
Curve 1 is a gaussian distribution for 
which a = .0028 radians. The distribu- 
tions are normalized to the same total 
number of events. From eq. (4), Chapter 
II/ the expected deflection angle is 
.0045 radians. 




PROJECTED ANGULAR DISTRiBUTION/MOMENTUM 
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MOMENTUM (MeV/c) x SCATTERING ANGLE (RADIANS) 

Fig. 15 . Calculated electron scattering-angle distribu- 
tion. Curve 1 is the angular distribution of 
the projected scattering angle. Plotted as 
above, the curve is valid for all electron 
momenta above a few MeV/c. Curve 2 is a 
gaussian distribution with a = .42 MeV/c 
radians. Only one sign of the scattering 
angle is shown since the distribution is 
symmetric about zero. 
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valid for all electron momenta of interest. Since p = R f 
the electron rigidity, we can substitute for p its value 
from eq. (4), Chapter II, to obtain 

a<K = — k_ e (6) 

3.55 

i.e., the standard deviation of the projected scattering- 

angle distribution is directly proportional to the "true," 

or magnetic, deflection angle e. Curve 2 -m Figure 15 is 

a gaussian distribution with a standard deviation of 

.42 MeV radians. At p<f> = .42 radians the integral 
c c 

distributions of curves 1 ,and 2 are equal. For calculating 
the deflection resolution P of the detector we shall use 
this gaussian curve as an approximation to the projected 
scattering-angle distribution. Since P for the distribu- 
tion is defined with the full width at half maximum, the 
true scattering-angle distribution would give a deceptively 
small value due to its high, narrow maximum at small angles. 
{Fifty percent of the particles have scattering angles 
greater 'than the value at half maximum for the true 
scattering-angle distribution; with the gaussian approxi- 
mation only 30 percent have larger scattering angles.) The 
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long large-angle scattering taxi will be taken into account 
in the trajectory consistency checking (see section B) . 

We can now derive the deflection resolution by com- 
bining the independent error distributions due to the 
intrinsic angular resolution and to multiple scattering. 
Using the gaussian approximation for the latter, we obtain 
the standard deviation erg of the deflection measurement 
for electrons by directly combining eqs. (4) and (6) to get 

Q 0 = (j, 2 + a A T ~ /(. 120) 2 + ( . 003) 2 , (7) 

We have used here k = .42 radians and a value for o A 

c A 

appropriate for both "perfect" and "one-error" events. The 
deflection resolution P, FWHM, is then 


2.36 On 

p ® V » 


(.18) 2 + (.0045/0) 


Equation (8) also applies for cosmic-ray nuclei, since 0 is 
very small for nuclei above the Cerenkov threshold and the 
scattering term (the first term under the radical sign) is 
therefore negligible. In Figure 16 we show cr A /0, cr^/®/ 
Oq/Q, and the resolution P. The measured values (triangles) 
for the resolution were obtained in a mono-energetic posi- 
tron beam at the California Institute of Technology 
Synchrotron. The deflection resolution P, as defined 



RIGIDITY (MV) 


Fig. 16. deflection resolution, P, FWHM, of the 
detector Data points (triangles) were 
measured in a mono-energetic positron ■ 
beam. Curves 1, 2 , and 3 are a^/0, 
Oi/0, and oq/0, respectively. 
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above , is derived for the distribution of measured deflec- 
tions about the "true," or magnetic, deflection angle 
which corresponds to the rigidity of the particle according 
to eq. (4), Chapter II. Thus P is a direct measure of the 
accuracy to which the detector can determine the particle 
rigidity. 

B) Event selection and detection efficiency 

Because all anti-coincidence counters are active and 
no interactions are involved in particle identification, 
events accepted for analysis are selected purely on the 
basis of good chamber operation and a self-consistent 
trajectory determination. The event analysis can be divi- 
ded into 3 steps: 1. determination of the trajectory in 

each spark chamber; 2. calculation of the deflection angle; 
and 3, checking whether a particle of the derived rigidity 
could travel a path consistent with the trajectory deter- 
mined by the spark 'chambers. 

Step (1): Initially each spark chamber box is 
examined independently. In each chamber a -least-squares 
fit of the measured spark locations is made to a straight 
line. Only those events are accepted where 3 or 4 of the 
gaps in each spark chamber show sparks within 2.5 mm of 
the best-fit straight line. Should one or more gaps of a 
4-gap fit not fulfill this criterion, the spark farthest 
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from the fitted trajectory is eliminated and a 3-gap fit 
is made. The average deviation of all gaps included in the 
final fit may not exceed 1.25 mm or the event is rejected. 
The mean deviation is generally less than 0.3 mm but tends 
to be slightly greater for trajectories away from the 
vertical. It is also required that the trajectory fit be 
within the acceptance cone of the detector. 

The original intention was to ignore those gaps 
where the multiple spark indicator (MSI) bit was set (see 
Chapter II, section B) . During flight several spark planes 
developed persistent spurious sparks at the edge away from 
the pick-up coil. Although these edge sparks did not seem 
to prevent registration of the trajectory spark location, 
they set the MSI, and thus made it necessary to ignore 
this bit. The effect is to increase the potential chance 
of fitting an incorrect trajectory in some ambiguous events 
or of failing to recognize a multiple particle event. The 
trajectory consistency checks, described in step (3) below, 
however, greatly decrease the probability of accepting an 
incorrect trajectory. Furthermore, a careful examination 

i 

by hand for possible incorrect trajectory assignments or 
multiple particle events showed very few questionable 
events (much less than one percent) . We therefore feel 
that ignoring the MSI bit has not significantly affected 


our results. 
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Step (2) : The deflection is calculated according to 

eq. (2), Chapter II, and is simply the difference of the 
entrance and exit angles. The particle is assigned a 
rigidity according to eq. (4), Chapter II. 

The alignment for zero deflection is determined by 
means of the cosmic-ray protons and alpha particles, which 
greatly outnumber other particle species triggering the 
detector and, in addition, are confined to small deflection 
angles. In Figure 17 we show the distribution of measured 
deflection angles 0, for |0|<_ .010 radians, measured 
during the nighttime period of flight C3. The criteria of 
step (3) below have been applied to these data. We have 
estimated the rigidity spectrum of cosmic-ray nuclei above 
the detector threshold by adjusting the observed spectra of 
Ormes and Webber' (57) (58) according to the Mt. Washington 

neutron monitor count rate measured during the flight. 
Curves 1 and 2 in Figure 17 are obtained by folding the 
estimated spectrum with an angular resolution for which 
Oq = .0027 radians and Og = .003 radians, appropriate for 
"perfect" and "one-error" events, respectively (eq. (4) . 

We conclude from the good agreement of the measured and 
calculated distributions that our deflection zero is accu- 
rate to within better than .0005 radians. Distributions 
measured during other flights are similar. 
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MEASURED DEFLECTION (MILLIRADIANS) 

Fig. 17. Distribution of measured deflection 0 with 
f 0 [ < .010 radians during nighttime period 
of flight C3. "Perfect" events (solid histo- 
gram) and "one-error" events (dashed histogram) 
are plotted separately. Curves 1 and 2 give 
the predicted distribution for cosmic-ray 
nuclei folded with angular resolutions for 
which 0 q = .0027 and c?g = .003, respectively. 
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Step (3) : The last step is to check the self- 

consistency of the trajectory, i.e., to determine whether 
a particle of the assigned rigidity and measured entrance 
trajectory could be expected to emerge from the magnet with 
the measured exit trajectory. An obvious procedure would 
be to mathematically propagate such a particle through the 
measured magnetic field in the same manner used in the cal- 
culation of the geometrical factor (Chapter II). If the 
mathematical particle emerged from the field with the 
measured exit parameters, we would accept the event. The 
relatively large scattering which can occur and which, of 
course, is random in nature, considerably complicates this 
approach , however . 

We use a much less time-consuming procedure, which 
utilizes the symmetry of the detector and the magnetic 
field and involves the detector resolution m a very 
straightforward manner. The trajectory parameters used in 
the checking are illustrated in Figure 18. We show there 
a projected view of a particle trajectory, assuming no 
scattering and an idealized magnetic field B which is uni- 
form and constant within the magnet gap and zero beyond. 

For such a field the path is an arc of a circle within the 

i 

field region and is straight above and below the magnet. 

I 

We construct the straight line ab between the trajectory 
positions at the bottom high voltage plane of the top 
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particle TRAJECTORY 



Pig. 18. Schematic view of a particle trajectory 
seen in projection, with definition of 
parameters used in trajectory self- 
consistency checking. 
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spark chamber and the top high voltage plane of the bottom 
spark chamber. Similarly, we construct the straight line 
cd between the trajectory positions at the top and bottom 
edge of the magnet gap. We define the angles X 2 , ^l • 

^2 and to as shown in Figure 18. Following the convention 
adopted in Chapter II for the angles a^, a 2 r and 0, clock- 
wise angles are taken to be positive. The following 
relationships can be derived directly from the figure and 
simple geometry. 


= - to - 


(9) 


X 2 = w + c*2 


( 10 ) 


0 a 2 a i “ X i + X 2 ~ + X 2 


x ' 2 - = 0 


( 11 ) 

( 12 ) 


Equation (12) is both a necessary and sufficient 

condition that an arc of a circle can be drawn between c 

i 

and d which will join smoothly to the trajectory external 
to the magnet gap. With the given field such a circular 

* f * 

arc is indeed an acceptable path for a particle of some 
rigidity. The presence of the non-uniform fringing field 
in the actual d.etector causes a real trajectory to deviate 
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from a circular arc, and, in general, also destroys the 
equality of A^ and A^ . However, we can define a new 
parameter 


A = 



(13) 


which will be useful for the trajectory checking. Deter- 
mination of A^ and- A^ only requires knowledge of the 
trajectory in the spark chambers. Using the measured' 
magnetic field, we have calculated the range of A for all 
possible valid trajectories (without scattering) within 
our acceptance cone and find that in all cases 


< .05 0 


(14) 


For most trajectories j A } is very much less than this 
limit. We shall see below that the measurement of A is 
subject to essentially the same errors as the measurement 

i 

of 0, which was discussed in section A. To the degree 
1 

\ > 

that A ,for the unscattered trajectory approximates zero, 

A is, in fact, a direct measure of the deviation of the 
measured entrance and exit path segments from a valid- 
unscattered particle trajectory through the detector. 

Combining eqs. (9) , (10) , and (13) .we have for A 


A = af + a 2 + 2w 


(15) 
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The standard deviation a ^ for the determination of A is 
therefore 



where cr u is the standard deviation of the measurement of 
a). The last step m eq. (16) follows from eq. (11). Since 
we have already included scattering in Oq, we need only 
consider the intrinsic angular error in a w . Points a and 
b of the trajectory are each defined by a pair of spark 
gaps. The derivation of a is essentially the same as that 
of a a (equations (1), (2), and (3)), with r / replacing r. 
There are two differences: 1. the selection criteria 

allow as few as 2 sparks to define the line ab (although 
in practice this occurs only rarely), and 2. the averaging 
over incident angles does not apply since w is restricted 
to a narrower range of angles than a-^ or a 2 . The standard 
deviation is given by 



.0011 radians for 4 good sparks 

(17) 

.0014 radians for 3 good sparks 
.0016 radians for 2 good sparks 
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since ct^ - .008" (eq. (5), Chapter II) and r 7 = 7". For 
our purposes it is sufficient to take an average value 

cx^ - .0014 radians (18) 

Substituting eqs. (7) and (18) into eq. (16), we obtain 

a A “ ^ ( . 120) 1 + (.003) 2 + 4 (. 0014) 2 

(19) 

~ ^ (.120) 2 + (.004) 2 

The content of eq. (19), for a is the same as that of 
eq. (7) for Oq , i".e., a term due to scattering and a term 
due to the intrinsic angular resolution of the instrument 
(though the latter is somewhat broadened) . Hence, in 
particular, the distribution of scattering angle shown in 
Figure 15 applies equally well to the A distribution. 

(Note, however, that the gaussian approximation is still 
used in eq. (19).) 'Comparing eqs. (14) and (19) we see 
that, within the resolving ability of the detector, A for 
all unscattered trajectories is very close to zero. Thus 
the measured value of A is a direct indication of the error 
in the trajectory determination. Furthermore, the distri- 
bution for A contains the angular resolution of the 
instrument in a straightforward way which makes it a good 
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parameter to determine whether given measured entrance and 
exit trajectories are segments of a self-consistent path. 

Based on these considerations and taking into 
account the true scattering-angle distribution, we choose 
as a criterion on a for acceptance of an event 

A 2 £ (.320) 2 + (.008) 2 (20) 

This choice of limits on A is sufficiently stringent that 
the probability is very small that an event will be 
accepted where the assigned trajectory includes spurious 
spark locations. On the other hand, only about 7 percent 
of all valid trajectories are rejected, independent of 
rigidity. This result can be seen from the integral 
scattering-angle distribution in Figure 15 (<f> = .320 
corresponds to p<f> = 1.15 in the figure) and from the 
gaussian error distribution with a = '.004 which applies for 
the intrinsic angular resolution. In applying this 
criterion we must use the measured 0, of course, while 
development of eq. (19) implicitly assumed the "true" 0, 
corresponding to the true rigidity of the particle. Because 
the resolution is quite .good for rigidities where the value 
of 0 is large, however, the dif f erehce-.does not produce a 
significant error. The criterion is valid for nuclei as 
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well as electrons since the electron scattering term is 
negligible for rigidities where nuclei trigger the 
instrument. 

Figure 19 shows A distributions for mono-energetic 
positrons of 800 and 100 MeV measured at the California 
Institute of Technology Synchrotron. The appropriate 
theoretical distributions, including the large-angle 
scattering, is also shown in each case. The rejection 
zones according to eq. (20) are indicated by cross- 
hatching . 

Essentially no bias according to particle species 
or rigidity has been introduced by the event acceptance 
criteria described in steps (1) and (3) above. Furthermore, 
the triggering requirements are such that all events are 
potentially acceptable, depending only on spark chamber 
performance. It follows that the detection efficiency 
D, valid for all rigidities, is simply 


D = 


n 

N 


(21) 


where N is the total number of coincidences and n is the 
number of accepted events. D was typically .7 - .8 
throughout all three flights. 
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Fig. 19. Measured distributions of the trajectory 
parameter A . 

Graph A. Mono-energetic beam of 800 MeV 
positrons. Mean deflection is 
.0045 radians (Compare Figure 
14) . 

Graph B. Mono-energetic beam of 100 MeV 
positrons. Mean deflection is 
.036 radians. 

Also shown in each graph is the theoretical 
distribution. The shaded areas are the 
rejection zones according to eq. (20) . 
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C) Dead Time Correction 

The phase A period is a natural time interval for 
evaluating the event rate. During such a 15-mmute period 
the total sensitive time t_ (m seconds) is given by 

t s = (900 - n e t e ) (1 - n a t a ) (22) 

where 

n e = number of events recorded 
t e = time required to write a word 
= .35 seconds 

n = total guard counter singles rate 

a (per second) 

t a = dead time following anti-coincidence 

puise = 2 x 10“ 6 seconds 

The above formula is correct provided n a is much 

5 

less than 5 x 10 counts per second, which was true through- 
out all flights. The value of used for a given period 
was the average of the rate recorded during the phase B 
periods immediately before and after the phase A period 
in question. 

At float altitude the fractional dead time was, 
typically 

t - t = 
s 


.17 t 
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The coincidence rate calculated using the sensitive time 
from equation (22) differs from the rate measured during 
the phase B period by less than one percent, which is well 
within the statistical accuracy of the phase B counter. 

D) Instrumental background 

The probability of a random pulse in the Cerenkov 

counter accompanying the passage of a particle below the 
v 

Cerenkov threshold is about 0.1 percent (see Chapter II, 
section A) . For heavy particles with rigidity less than 
200 MV this represents a negligible contribution to the 
measured flux. 

Possible sources of contamination from particles 

1 

v 

above the Cerenkov threshold are discussed individually 
below. 

1. Nucleonic component 

Cosmic-ray protons and alpha particles represent 
about 95 percent of all particles triggering the detector 
system at float altitude. Because they are effectively 
confined to small deflection angles, however, their contri- 
bution to rigidities below 200 MV can only result from 
secondary particles or from a possible miscalculation of 
the trajectory. The probability of the latter was 
measured in a high energy (1000 MeV) electron beam. The 
data showed fewer than 0.5 percent of all events which 
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appeared to have rigidities less than 200 MV. A similar 
percentage for nuclei in flight would give a rate of 
3 x 10 counts per second^ this is only about 5 percent 
of the rate in any one of our 5 rigidity intervals. This 
limit includes any contribution due to knock-on electrons 
produced in the detector. The total material of the 
detector in the beam above or within the spark chambers 
is about .005 interactions lengths; hence contamination by 
particles resulting from nuclear interactions in the 
detector is also not significant. 

2. Pions and muons 

Pions and muons with momentum greater than about 

140 MeV/c and 105 MeV/c, respectively, are above the 

v 

threshold of the lucite Cerenkov counter. At our float 
altitude pions near 200 MeV/c decay within .002 g/cm 2 of 
their point of production (59) . Hence their contribution 
relative to muons may be ignored. The muon spectrum at a 
depth x << L a is, from Verma (59) , 


F(E y ,x)a Ey « F t (E t ) 


as 


IT 

as 

y 



X 

L 


f 

_1 x 1_ 

g+1 L a g+2 

» j 


(23) 


1 * 

Here P 7r (E^) is the charged pion production spectrum, 
and Ep the pron and muon energies, L x and the proton 
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interaction and absorption lengths, respectively, and 
g = 1 . 07 /Ejj , where is the muon energy in BeV. We 
have calculated the muon spectrum at an atmospheric depth 
of 2.4 g/cm , using the pion production spectrum of Perola 
and Scar si (60) corrected to the proton spectrum at Fort 
Churchill m 1968. We make the usual assumption that the 
pion and resultant decay muon move at the same velocity 
and use L a = 120 g/cm 2 and = 100 g/cm 2 . The calculated 
flux of muons with rigidities between 100 and 200 MV is 
~10 -2 muons/ (m 2 * sec* sr. MV) . This is about 5 percent of 
the total measured flux in this interval. 

i 

3 . Gamma-rays 

v 

Gamma-rays may interact in the lucite Cerenkov 
radiator, trigger that counter, and subsequently scatter 
upward into the acceptance cone of the detector (see Figure 

t 

20) . In order to determine the sensitivity to this effect 
we have calibrated the detector in a y-ray beam at the 
California Institute of Technology Synchrotron. Measure- 
ments of the high-altitude y-ray spectrum at various zenith 
angles have been made by Fichtel, Kniffen and Ogelman (61) . 
Their observations were made at Mildura, Victoria, 
Australia, m December, 1966. The geomagnetic cutoff 
rigidity is 4.7 GV at Mildura (62), which is much higher 
than the value at Fort Churchill (<200MV) . The level of 
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e + e 



Fig. 20. Schematic representation of v gamma-ray 
interactions in the lucite Cerenkov 
counter which can produce spurious elec- 
tron events. The interactions involved 
are Compton scattering and pair produc- 
tion. Path marked e + will be registered 
as e“ and vice versa due to the backward 
trajectory. 
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solar modulation m 1966 was considerably less than in 
1968, however. The pion production curves of Perola and 
Scar si (60) indicate that these two effects are roughly 
compensatory. We have therefore used the y-ray spectra of 
Fichtel, et a_l. for Fort Churchill m 1968. Figure 21 
shows the derived contribution from this source to our 
measured positron and negatron spectra at float altitude 
(63) . 

With the exception of atmospheric secondary elec- 
trons, the y-ray produced background is the most 
significant contamination of our measured intensities at 
low energies. The contribution to the positron flux 
between 6 and 12 MeV is about 20 percent of the derived 
primary positron intensity. For negatrons and at higher 
energies the relative contribution is considerably less. 

The large error limits which are shown in Figure 21 are 
estimated errors based on the quoted accuracy of the y-ray 
spectrum at balloon altitudes and on the statistical and 

experimental limitations of the machine calibrations. 

\ 

Further calibrations are being .undertaken to improve our 
knowledge of the sensitivity of our detector to y-ray 
produced background. The corrections do not materially 

r 

affect our final conclusions regarding the primary electron 


spectrum, however. 



DIFFERENTIAL FLUX (M 2 -SEC-SR-MeV) H 



KINETIC ENERGY (MeV) ' KINETIC ENERGY (MeV) 

Fig. 21. Gamma-ray produced contribution to measured positron and negatron 
differential kinetic-energy spectra at balloon altitudes. 
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4. Splash albedo electrons 

In Figure 22 we show the splash albedo spectrum 
measured at Fort Churchill in 1967 by Israel (18) . These 
particles result primarily from very high energy primary 
cosmic rays and hence the spectrum m 1968 should be 
similar. Curve 1 shows this spectrum after propagation 
through the -15 g/cm^ of material in the lower part of 
our instrument. The lucite Cerenkov counter has an 
efficiency of . 04 for the detection of backward-moving 
particles. Multiplying curve 1 by the backward detection 
efficiency, we obtain curve 2 for the contribution to our 
measured electron spectrum from this source . The splash 
albedo should consist of essentially equal numbers of 
positrons and negatrons. The contribution to our measured 
flux is less than the statistical accuracy throughout our 
energy interval. 

E) Atmospheric Secondaries 

A significant part of the electrons between 6 and 

2 . 

200 MeV observed at 2.4 g/cm is of atmospheric origin. 

The atmospheric electrons result primarily from the decay 
of pions produced in interactions of the cosmic-ray nuclei 
with air nuclei; below 20 MeV knock-on electrons also con- 
tribute significantly. At somewhat greater atmospheric 
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Fig. 22. Differential kinetic-energy spectrum of 

splash albedo electrons. The measurements 
of Israel (Ref. 18) were made m 1967. 

Curve 1: the spectrum after propagating 

through 15 g/cm^ of material m 
the lower gondola. 

Curve 2 : the background contribution to 

our measured spectrum of down- 
ward-moving electrons, obtained 
by multiplying curve 1 by the 
backward detection efficiency. 


. 04 . 
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Figure 22 
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depths electromagnetic cascade showers become the 
dominant source. 

Calculations of the atmospheric secondary electron 
spectrum have been published by Verma (59) and by Perola 
and Scarsi (60) . The two calculated spectra differ 
considerably below several hundred MeV. In addition, 
neither includes knock-on electrons nor attempts to separ- 
ate positrons and negatrons. A new calculation has been 
made recently by Beuermann (47) . This calculation uses 
the pion production spectrum of Perola and Scarsi (60) 
corrected to the incident proton spectrum m 1968; it also 
includes knock-on electrons. The spectra of atmospheric 
secondary positrons and negatrons are calculated separa- 
tely over a wide range of atmospheric depths. The depth 
dependence of the residual primaries can be calculated 
in a similar manner, given an incident spectrum. 

We separate the primary and secondary contributions 
to our measured flux in the following manner. The atmos- 
pheric depth dependence of the total positron or negatron 
flux J * (d) , for the i'*'* 1 energy interval, is of the form 

J i ± (d) = a i ± s i ± (d) + b i ± p i ± (d) (23) 

where d is the -atmospheric depth, s i ± (d) and pj_“ (d) give 
the functional form of the calculated depth dependence of 



- 86 - 


the flux of secondary and primary positrons or negatrons, 
respectively, and and b^~ are parameters giving the 

relative contribution of each component. We determine the 
parameters a^ and by making a least-squares fit to 
seven data points from 2.4 to 42 g/cm 2 atmospheric depth. 
Since p^- (d) depends on the unknown incident primary spec- 
trum, an iterative process must be used. The derived 
primary contribution at float altitude, 

t i o 4 

b^~p^“ (d = 2.4 g/cm ) , is not very sensitive to the choice 
of the incident primary spectrum, p ± (d = 0), however, and 
so the iteration converges quickly. 

In Figure 23 we show the measured growth curves and 
the residual primary and atmospheric secondary contribu- 
tions determined by the least-squares fit. Included are 
data from the ascent (which in each case occurred after 
the evening transition to lower geomagnetic cutoff) and 

from the nighttime interval of the float period. (See 

/ 

Chapter V, section A for a discussion of the diurnal 
cutoff variation) . The data from the 3 flights showed no 
systematic variations and have therefore been combined for 
greater statistical accuracy. The energy intervals 
indicated are those at the top of the detector. The x z 
probability P for each fit is included m the graphs. In 
the 50-100 MeV interval, the calculated atmospheric 
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secondary growth curves for both positrons and negatrons 
fit the measured curves within statistical errors. The 
fit for both results in a small negative primary contri- 
bution. For this interval, therefore, only upper limits 
can be derived. The measured spectra of positrons and 
negatrons together with the separation into primary and 
secondary components resulting from our fitting technique 
are shown in Figure 28 and discussed in Chapter V, section 
B. 
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Fig. 23. Measured positron and negatron event 
rates vs. atmospheric depth for the 
nighttime period. The energy intervals 
indicated are those measured at the 
top of the detector. Also shown is 
the separation into primary and 
atmospheric secondary components as 
determined by the least-squares 
fitting technique described m the 
text. 

Dashed curve: best-fit primary contri- 

bution 

Dash-dot curve: best-fit secondary 

contribution 

Solid curve: best-fit total positrons 

or negatrons 

The x 2 probability, P, is indicated for 
each fit. The increase in the intensity 
of the primary component with increasing 
depth, which is evident in some of the 
graphs, is due to the assumed form of 
the extraterrestrial positron and 
negatron spectra (see Fig. 28) . 
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V. RESULTS AND DISCUSSION 


A) Diurnal variation of the geomagnetic cutoff 

Fort Churchill is the northernmost balloon- launching 
station where the facilities and terrain permit experiments 
to be handled with relative ease and good probability of 
equipment recovery after flight. The nominal geomagnetic 
cutoff rigidity, calculated for the internal geomagnetic 
field only, is ~200 MV at Fort Churchill (62) . The obser- 
vation of primary cosmic rays considerably below this 
rigidity is made possible by a large decrease in the local 
geomagnetic cutoff rigidity at night. This diurnal cutoff 
variation is of great interest in itself as a tool to 
study the magnetosphere. Our concern here, however, will 
be with the variation as it relates to low-rigidity primary 
cosmic rays. More general discussions can be found in 
the literature. 

1 . Background 

The geomagnetic field prevents low-rigidity cosmic 
rays from reaching the earth at low latitudes. Calculations 
of the geomagnetic cutoff rigidity based on the internal 
field of the earth have been performed by several investi- 
gators. At the same time it has become apparent that the 
actual cutoffs differ from these calculations, particularly 
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at high latitudes. In addition, a diurnal variation in 
the local cutoff rigidity has been observed at geomagnetic 
latitudes between 65° and 75°. This cutoff variation has 
been observed in polar cap absorption (PCA) events (64), 
in measurements of low energy solar proton fluxes made on 
polar orbiting satellites (65) (66) , and in measurements 

of low energy electrons observed with balloon-borne 
detectors (11) (46) (67) . The electron data have shown 

that the cutoff near Fort Churchill (invariant latitude 
=70°) varies from above 100 MV during the day to below 
10 MV at night. The transitions, which last 1-2 hours, 
occur at about 0600 and 1800 local time. 

A consistent theoretical picture has emerged which 

i 

attributes the cutoff variation to the influence of the 
geomagnetic tail (68-71) . In Figure 24 we show schemati- 
cally a model of the magnetosphere constructed by Williams 
and Mead (72) based on available satellite data. At 
geomagnetic latitudes £ 68° the field lines retain a quasi- 
dipole shape at all times, while at high latitudes, >_ 80°, 
the field lines are always swept back into the tail. At 
intermediate latitudes, 68°-80°, however, the field lines 
change from a near dipole shape during the day to extension 

t 

deep into the tail at night. There is evidence that 
interplanetary particles have essentially free access to 



Pxg. 24. Geomagnetic field lines in the Williams and Mead 
model magnetosphere. The horizontal scale 
indicates distance in units of earth radii. Indi 
vidual field lines are labeled by the geomagnetic 
latitude at which they intersect the earth. 
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the tail, where they can travel freely along the field 
until they reach the earth (73-75) . The cutoff at the 
foot of a geomagnetic field line which extends into the 
tail should therefore be essentially zero. On the other 
hand, if the field line is a quasi-dipole line, connecting 
directly with a conjugate point in the opposite hemisphere, 
one would expect the local cutoff to be similar to that 
calculated using only the internal field. Detailed calcu- 
lations of particle orbits in the Williams and Mead model 
magnetosphere have been performed which support the general 
picture outlined above (70-72) . 

This theoretical model leads to the following 
interpretation of the electron measurements (46) . During 
the day the cutoff at Fort Churchill is somewhere between 
about 100 and 200 MV. The electron flux measured below the 
cutoff is a mixture of atmospheric secondaries and return 
albedo. At about 1800 local time the local field lines 
are swept back into the tail and the cutoff decreases to 
near zero. The return albedo particles are then replaced 
by the interplanetary electron flux, which is of lower 
intensity. The measured flux therefore decreases sharply. 
At about 0600 local time the reverse transition takes 
place. This interpretation of the measured electron flux 
variation is supported by the data of Israel (18) who 
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measured the splash albedo spectrum at Fort Churchill and 
found it to be in good agreement with his measured 
downward-moving daytime electron spectrum below about 
100 MeV. 

2. Results and discussion 

We have referred previously to the timing of our 
observations which enables us to use data gathered during 
ascent to directly correct our nighttime measurements for 
the contribution of atmospheric secondary electrons (Chap- 
ter IV, section E) . We present here profiles of count 
rate versus local time which show the location of the 
morning transition during each flight and the limits of the 
nighttime and daytime periods' used in subsequent analysis. 
We also present evidence that the nighttime cutoff was 
below our analysis threshold of 6 MeV at the detector. In 
addition, our measurement of the charge ratio of the day- 
time flux is entirely new evidence supporting the albedo 
origin of this component. 

In Figure 25, we show our hourly count rate plotted 
against local time for each of the three flights. Positron 
(dotted histogram) and negatron (solid histogram) count 
rates are shown separately for 5 energy intervals between 
6 and 200 MeV at the detector. Data from 4 consecutive 
phase A periods, an elapsed time of 64 minutes, are 
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Fig. 25. Event rate vs. local time. Negatrons (solid 
histogram) and positrons (dotted histogram) 
are shown separately. Typical 1-a error limits 
are indicated. The limits of the time intervals 
included in summaries of "daytime" and "night- 
time" data are indicated. Also shown are the 
contributions from atmospheric secondary 
negatrons (long dashes) and positrons (short 
dashes) . In the energy intervals 25-50 MeV 
and 50-100 MeV positron and negatron secon- 
dary intensities are approximately equal; 
therefore only one component is shown. Varia- 
tions in the secondary rate reflect altitude 
variations . 

Graph A. Flight Cl 
Graph B. Flight C2 
Graph C. Flight C3 

A period of bad data near 0800 local time 
during flight C3 has been omitted. 
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averaged for each point. Typical error limits are indi- 
cated. The dashed curves (long dashes for negatrons and 
short dashes for positrons) give the contribution of 
atmospheric secondaries as determined by our least-squares 
fit. 

In flight C2 the morning transition occurs later 
than for flights Cl and C3 (about 0700 local time versus 
approximately 0300 and 0430 local time) . It is also more 
abrupt. The trajectories of all three flights were simi- 
lar (Figure 11) so that the difference cannot be due 
simply to a different location of the detector in the geo- 
magnetic field. A small solar event occurred during the 
latter part of flight C2, around 1600 UT (0800 local time 
at the detector) , but it is difficult to imagine a connec- 
tion with the delayed transition which actually preceded 
the solar event. Large variations in the time and abrupt- 
ness of the transitions have also been observed by Israel 
and Vogt (67) m 1967, including one flight (their flight 
Cl) where the morning transition had a time profile similar 
to the transition seen in our flight C2. On the other hand, 
they did not see a transition as early as the one during 
our flight Cl in any of their three flights. All of the 
flights of Israel and Vogt occurred during quiet geomag- 
netic conditions. Jokipii, L'Heureux, and Meyer (44) 
observed 5 transitions, all between 0500 and 0700 local 

i 
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time. Hovestadt and Meyer (76), flying within several 
weeks of our flights, observed a transition at -0400 local 
time. On the other hand. Rocks tr oh and Webber (12) , whose 
flight was only 2 days after our flight C3, claim to see 
a morning transition after -0900 local time (77) . It is 
possible that day-to-day changes in the configuration of 
the magnetosphere are causing the large variations in the 
transition time. Despite the fact that there is no 
statistically significant change in the measured electron 
flux in flight C2 until the transition near 0700 local 
time, we have included in our nighttime data interval only 
the period prior to 0945 UT (0304 local time) . We thus 
make the nighttime period included in the data summary 
similar for all flights. 

The energy interval between 100 and 200 MeV displays 
a more diffuse transition than the other intervals in all 
flights. During flight the balloon moved westward, passing 
to lower geomagnetic latitudes (see Figure 11) , and conse- 
quently to higher dipole cutoffs. We are thus probably 
seeing a mixture of the initial transition from the tail 
field cutoff near zero MV to a closed field cutoff near or 
somewhat above 100 MV followed by a more gradual change 
due at least partly to the motion of the balloon. 

In Figure 26 we show the measured daytime and night- 

i 

time positron and negatron raw spectra for each flight. 
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Fig. 26. Measured differential kinetic-energy spectra 

of positrons and negatrons for daytime and 

nighttime intervals. The data from each flight 

are plotted separately. Nighttime spectra are 

' 2 

for an average floating depth of 2.4 g/cm 
residual atmosphere. Daytime spectra are for 
an average depth of 2.2 g/cm . The daytime 
data of flight C3 have been corrected for the 
change in altitude toward the end of that 
flight. All other data are raw fluxes. 


Graph A. 
Graph B. 
Graph C. 
Graph D. 


Positrons - night 
Negatrons - night 
Positrons - day 
Negatrons - day 
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Figure 26 B 
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(The daytime data for flight C3 have been corrected for 
the increased atmospheric secondary flux caused by the 
change in altitude toward the end of that flight. This is 
done in order that all daytime measurements apply to an 

9 

average depth of 2.2 g/cm residual atmosphere. All other 
values shown, day and night, are raw fluxes) . Within the 
statistical uncertainty, the agreement among flights is 
good . 

Comparison of the day and night spectra reveals 
large differences down to the lowest energy interval, 

6-12 MeV at the detector. Below 6 MeV the geometrical 
factor of the detector goes quickly to zero and the elec- 
tron detection efficiency of the Cerenkov counter is also 
declining. It is therefore not possible to use this 
interval for a reliable flux determination. Nevertheless, 
comparison of the measured nighttime and daytime rates, 
given m Table 4, shows that there is a well-defined 
diurnal variation for these electrons also. We conclude 
that the nighttime cutoff is well below our detector 
threshold of 6 MeV (corresponding to 11.1 MeV at the top 
of the atmosphere) , and that our nighttime data are there- 
fore free of contamination by return albedo. 

In Figure 27 we show our measured daytime positron 
and negatron spectra, as well ass their sum, corrected to 
the top of the atmosphere. We combine here the data from 
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TABLE 4 


MEASURED RATES OF POSITRONS AND NEGATRONS BELOW 6 MeV 

Rates for each flight and for daytime and 
nighttime intervals are given separately. 

The quoted errors include only statistical 
uncertainties . 


Nighttime Daytime 

Flight Component Rate (hr 1 ) Rate (hr -1 ) 


+ . 6 


+ 2.3 
- 1.5 


+ 1.0 

- . 6 


+ 2.3 

- 1.5 


+ . 9 
- .4 


+ 1.6 


+ 1.1 

- .7 


+ 2.3 
- 1.7 


+ .2 

- .3 


+ . 8 

- . 6 


+ 1.3 

- 1.0 


+ 1.1 

- 1.0 
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KINETIC ENERGY (MeV) 

Fig. 27. Daytime differential kinetic-energy spectra of 
positrons, negatrons, and total electrons, 
incident at the top of the atmosphere. The data 
of Israel (Ref. 1 (18)) are for the splash albedo 
in 1967. 
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all three flights. Also shown is the splash albedo spec- 
trum at Fort Churchill measured by Israel (18) in 1967. 
Since the splash albedo flux is produced primarily by the 
nucleonic component above several BeV we would expect the 
spectrum in 1968 to be almost the same as in 1967. Our 
daytime flux agrees well with the splash albedo of Israel. 
In addition, the similarity of the spectra of the positron 
and negatron components can be taken as new evidence of the 
albedo origin of the daytime electron flux. Equal contri- 
butions of positrons and negatrons are expected in the 
albedo spectrum, which should be similar to the spectrum of 
atmospheric secondary electrons deep m the atmosphere. On 
the other hand, the spectra of extraterrestrial positrons 
and negatrons measured during the nighttime interval differ 
considerably from one another (see the next section.) 

B) Primary cosmic-ray positrons and negatrons 
1 . Background 

Our measurements of the cosmic-ray positron and 
negatron spectra have bearing upon such astrophysical 
topics as solar modulation and the origin and interstellar 
spectrum of the electron component of the cosmic rays (see 
Chapter I). In this section we shall briefly review the 
current state of theoretical and experimental knowledge of 
solar modulation as it pertains to our new data. A general 
review can be found in Webber (78) . 
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While no theory of solar modulation has successfully 
explained all features of the observed temporal variations 
of the spectra of cosmic rays, there is reasonable agree- 
ment with the diffusion-convection model, first advanced 
by Parker in 1958 (79). In simple form, this model 
relates j E (R,t), the differential rigidity spectrum of a 
given cosmic-ray constituent observed at the earth at time 
t, to the interstellar spectrum j 0O (R) by 


j E (R/t) 


j«,(R) exp 



n(t) 

6 f(R) 


( 1 ) 


Here n is a time dependent but rigidity independent para- 
meter related to the geometry of the modulation region, 

8 is the particle velocity in units of the speed of light, 
and f (R) is a function of the particle rigidity R. 

According to the work of Jokipii (80) (81), f (R) can be 

related in a straightforward way to the power spectrum of 
the magnetic irregularities in the solar wind. Measurements 
of the power spectrum imply (82) 


f (R) 


! R" 5 for R > R 

° 

R for R < R 

o o 


( 2 ) 


where 6 is between 0.5 and 1. The transition rigidity R q 
is related to the correlation length of the magnetic 
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lrregularities in the solar wind and should be on the order 
of several hundred MV based on the measurements of the 
magnetic field (82) . 

The above formulation of the diffusion-convection 
theory neglects the effect of adiabatic deceleration during 
passage of the cosmic rays through the solar wind. The 
effect is included in the Fokker-Planck equation proposed 
by Parker (83), but a general solution is difficult, and 
up to this time only special cases have been examined in the 
literature (81) (83-86) . Gleeson and Axford (87) and Fisk 

and Axford (88) have used a somewhat different approach to 
the same problem. In all cases, however, the application 
at low rigidities is uncertain due to a lack of knowledge 
of the applicable diffusion coefficient. 

The most extensive cosmic-ray data exist for protons 
and helium nuclei. Because knowledge of their interstellar 
spectra is lacking, however, it is not possible to determine 
the absolute modulation, and so study of the solar modula- 
tion of protons and helium nuclei has been necessarily 
limited to temporal variations. From equation (1), it 
follows that the relationship of the intensities at two 
different times, t^ and t 2 , is 

b m = n(t 2 ) - nltx) 

: E ( R ' t 2> 


f (R) 


( 3 ) 
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provided there is no change in the functional form of f (R) . 
Time variations of the proton and alpha particle spectra 
give reasonable agreement with the functional form of f (R) 
given in eg. (2) with <5 ^ 0.5-1. 0 and R Q - 500 MV (78) (82). 

The data do not extend below several hundred MV rigidity, 
however . 

Several investigators have studied temporal changes 
in the measured electron spectrum and have come to differ- 
ing conclusions. L'Heureux, Meyer, Verma, and Vogt (89) 
report essentially no change m the electron spectrum be- 
tween 250 and 1050 MV in the period 1960-1966. They quote 
an upper limit of 60 percent for the fractional change in 
modulation during this time. Bleeker, Burger, Deerenberg, 
Scheepmaker, Swanenburg, and Tanaka (90) report no change 
above 500 MV between 1965 and 1966. On the other hand, 
Rockstroh and Webber (12) claim to see significant continu- 
ing modulation from 1965 through 1968 for rigidities 
between 20 MV and 1 GV. In addition, Bleeker, et al. (13) 
have recently reported a reduction in the electron intensity 
between 500 MV and 5 GV of about 30 percent between 1966 
and 1968. However, the actual reported absolute fluxes of 
Bleeker, et al. differ considerably from those of Rockstroh 


and Webber. 
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In view of these apparent inconsistencies, the 
question of the existence of a temporal change in the elec- 
tron spectrum must be considered unresolved at the present 
time. This question is an important one. The temporal 
variation of the proton and helium modulation is well* 
known. Electrons, with their much smaller mass-to-charge 
ratio, could be of great value in distinguishing between 
possible models of the modulation mechanism, as well as 
extending the measurements to very low rigidities where 
measurements of nuclear species are particularly difficult. 

The problem of ultimate concern, however, is the 
determination of the total, or absolute, solar modulation, 
for which it is necessary to know or be able to estimate 
the interstellar spectrum j^CR) as well as the directly 
measured spectrum at the earth, j (R,t). This restriction 
severely limits the cosmic-ray constituents which can be 
used to determine the absolute solar modulation. Ramaty 
and Lingenfelter (91) have calculated the interstellar 
spectra of deuterium and helium-3 assuming that these par- 
ticles are produced in collisions of cosmic-ray nuclei with 
the interstellar gas. Comparing with measured intensities 
at rigidities > 600 MV they found good agreement with the 
modulation function of eqs. (1) and (3) with 
n = 350 ± 150 MV, 5 - 1 , and R q < 600 MV. Experimental 
difficulties prevent, reliable measurements at lower 
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rigidities. Several authors have attempted to deduce the 
interstellar electron spectrum from the galactic background 
synchrotron radiation (92-94) . These calculations involve 
the estimation of certain galactic parameters, notably 
the interstellar magnetic field. In general, the results 
indicate modulation m reasonable agreement with eg. ( 1 ) 
and (3) with n = 500 - 1000 MV and R Q 4 300 MV. Inter- 
stellar absorption below - 10 Mhz limits this method to 
rigidities > 200 MV. 

This restriction to rigidities above several hundred 
MV does not apply to the derivation of the absolute solar 
modulation employing cosmic-ray positrons. We make the 
assumption that all of the cosmic-ray positrons with 
energies above a few MeV originate in the decay of charged 
pions produced m interstellar collisions between high- 
energy cosmic-ray nuclei and the ambient matter. The 
interstellar negatron and positron spectra from this 
source can be calculated with reasonable accuracy (27-29). 
We shall use here the calculation of Ramaty and Lingen- 
felter (28), which uses an integral path length of 4 g/cm^ 
for cosmic rays in the interstellar medium. By comparing 
our measured positron spectrum at the earth with the calcu- 
lated interstellar spectrum, we derive the absolute solar 
modulation of positrons between 11 and 200 MV. The 
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conditions under which this modulation can be applied to 
other cosmic-ray constituents, m particular to negatrons, 
will be examined later. i 

t 

2. Results 

t 

In Figures 26A and 26B we show the positron and nega- 
tron spectra measured during the nighttime interval of each 
of our three flights. In general the mutual agreement is 
good within the statistical accuracy of the measurements. 

In particular, there are no systematic deviations which 
might be associated with the recovery phase of the Forbush 
decrease which is evident in the neutron monitor count 
rates shown m Figure 13. The altitude profiles of all 
three flights were very similar to well past the time of 
the morning cutoff transition (Figure 12) . The average 
nighttime float altitudes were 2.45, 2.40 and 2.35 g/cm 2 
for flights Cl, C2, and C3, respectively, with variations 
of ± .15 g/cm 2 . We therefore consider the 3 flights 
equivalent and have combined the data for improved statis- 
tical accuracy in the subsequent analysis. 

In Figure 28 we show the combined measured nighttime 
spectra at 2.4 g/cm 2 residual atmosphere. We also show the 
division into primary and secondary components according 
to the least-squares fitting technique described in Chapter 
IV, section E. The fits are made to the raw data, i.e., no 
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Fig. 28. Nighttime differential kinetic-energy spectra of positrons and negatrons 
at 2.4 g/cm^ residual atmosphere. Shown are the measured spectra and 
the separation into primary and atmospheric secondary components as 
determined by the least-squares fitting technique described in the text. 
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corrections have been included for deflection resolution or 
for background events. The effect of folding in the resolu- 
tion of the instrument was calculated and was found to 
influence the measured flux by less than one percent m any 
energy interval. This correction is not included m the 
quoted data. The atmospheric depth dependence of the prin- 
cipal causes of background events - y-rays, splash albedo, 
and muons - is such that each source contributes almost 
exclusively to either the "primary" or the "atmospheric 
secondary" component as determined by the fit. In the 
case of the y-ray and splash albedo background the contri- 
bution is to the separated "primary" spectrum. Background 
subtraction after the fit has been adopted in order to 
facilitate possible later adjustments in the background 
corrections should new data become available. This is 
particularly important for the y-ray corrections where the 
present knowledge, both of the y-ray spectrum at float 
altitude and of our detector sensitivity, is somewhat 
limited. 

The "primaries" determined by the fit are 20-30 
percent or less of the total measured flux with the single 

I 

exception of the 6-12 MeV positron point where the contri- 
bution is ~ 60 percent. Statistical error limits are 
correspondingly large. All error limits shown in Figure 28 
are 1-a statistical errors only. The fitting procedure 
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determmed small negative "primary" contributions in the 
50-100 MeV interval for both positrons and negatrons (see 
Figure 23). Hence we can indicate only upper limits for 
this interval. The upper limits shown m Figure 28 for 
the 50-100 MeV interval are the 1-a uncertainties on the 
"primary" fit, taken from zero flux. 

After subtraction of the y-ray and splash albedo 
background the primary spectrum is corrected for energy 
loss to the top of the atmosphere. Only ionization loss is 
included. Bremsstrahlung loss also affects our highest 
energy interval. The magnitude of the correction depends on 
the spectral form at somewhat higher energies, however, and 
there is presently no general agreement on the spectrum 

I 

immediately above 200 MeV. In any case the effect on the 
primary flux between 100 and 200 MeV is less that ± 5 per- 
cent for any reasonable extension of our own spectrum to 
higher energies. This value is small compared with the 
statistical errors. Our derived primary positron and 
negatron spectra at the top of the atmosphere are shown in 
Figure 29. The steps leading to 'our incident spectra are 
summarized in Table 5. 

3. Discussion 

In Figure 30 we show our primary electron (e + + e ) 
spectrum together with some recent results of other 
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Fig. 29. Differential kinetic-energy spectra of positrons 
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TABLE 5 


NIGHTTIME POSITRON AND NEGATRON FLUXES 

_ 2 — 1 

A) Positrons. Flux values are in <m -sec’sr’MeV) 


Rigidity interval at 
the detector (MV) 

6-12 

12-25 

25-50 

50-100 

Flux from combined 
flights 

. 71± . 13 

.441.06 

.241.03 

.131.015 

Least-squares fit 
separation 

1) "Primaries" 

2) "Secondaries" 

. 44± . 13 
. 26±. 05 

.151.07 

.281.04 

.0301.041 
. 211.03 

-.0211.026 

.151.02 

Background corrections 
1) Splash albedo (1) 

. 03±. 01 

. 031.01 

.0161.005 

.00451.0045 

2) y-ray (1) 

. 09±. 09 

. 0441.030 

. 0061.003 

— 

Corrected primary 
positrons at 2.4 g/cm 2 
residual atmosphere 

.311.16 

. 0781.069 

.0071.042 

0+. 026 

Energy interval 

at top of atmosphere 

(MeV) 

11.1-17.1 

17.1-30.2 

30.2-55.2 

55.2-104.9 

Flux at top of 

.311.16 

. 0771.068 

. 0071.041 

0+.026 


atmosphere 


100-200 
. 0751.008 

. 010± . 012 H 

. 065± . 009 m 

i 

.00101.0003 
.009+. 012 
104.9-203 . 9 
. 0091.012 


(1) Error limits are estimated errors 






TABLE 5 (Cont.) 


B) Negatrons. 


Flux values are in 



sec • sr *MeV) 


1 


Rigidity interval at 
the detector (MV) 

6-12 

12-25 

25-50 

50-100 

100-200 

Flux from combined 
flights 

2. 80±0. 30 

.81+. 08 

.30+. 03 

. 11± . 02 

.084+. 008 

Least-squares fit 
separation 

1) "Primaries" 

. 72± . 40 

. 2 6± . 10 

. 072± . 047 

-. 038±. 030 

.030+. 014 £ 

2) "Secondaries" 

2 . 05± . 29 

. 55± . 07 

. 22± . 03 

.15±. 02 

. 051± . 009 f 

Background corrections 
1) Splash albedo (1) 

.03+. 01 

.03+. 01 

.016+. 005 

.0045+. 0045 

.0010+. 0003 

2) y-ray (1) 

.10±. 06 

, 0 2 5 ± . 0 1 0 

. 0066± . 0030 

— 

— 

Corrected primary 
negatrons at 2.4 g/cm^ 
residual atmosphere 

.59+. 41 

. 20± .10 

. 05 0± . 048 

0+.030 

.0291.014 

Energy interval at 
top of atmosphere 
(MeV) 

11.1-17.1 

17.1-30.2 

30.2-55.2 

55.2-104.9 

104.9-203.9 

Flux at top of 
atmosphere 

.58+. 40 

.20+. 10 

.049+. 047 

0+. 030 

.0291.13 


(1) Error limits are estimated errors 
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KINETIC ENERGY (MeV) 

Fig. 30. Differential kinetic-energy spectrum of 

extraterrestrial cosmic-ray electrons (e + + e“) 
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mvestigators . The data of Simnett and McDonald (21) and 
of Fan, L'Heureux, and Meyer (17) are satellite measure- 
ments and are thus free of contamination by atmospheric 
secondaries. Where their measurements overlap our own 
the agreement is very good. This is a good indication of 
the validity of our correction technique for atmospheric 
secondaries . 

Earlier measurements of the electron spectrum in 
the energy range 10-200 MeV are generally higher than our 
values measured in 1968; this is true even if restricted 
to satellite measurements and to balloon data where the 
local cutoff is known to be below the detector threshold 
throughout the period included in the analysis (for a 
summary of such measurements see Israel (18)). Some of 
the variation might be due to increased solar modulation 
m 1968, but differing atmospheric secondary corrections 
may also play a significant role. The large contribution 
of atmospheric secondaries which are present in balloon 
data covering this energy interval render the derived pri- 
mary spectrum very sensitive to this correction. The 
reduction technique for earlier satellite measurements at 
MeV energies has also been reexamined (21) . The deriva- 
tion of temporal changes of the solar modulation of elect- 
trons by combining data of different investigators is thus 
open to serious question and' we shall not attempt it here. 
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In Figure 31 we show the positron fraction of the 

total electron flux, N , / (N , + N -) as a function of 

e+ e+ e 

energy. Plotted are our values between 11 and 204 MeV 
together with those of Fanselow, Hartman, Hildebrand, and 
Meyer (35) at higher energies. A positron fraction for 
the interval between 55 and 105 MeV cannot be derived be- 
cause we have only upper limits for both components. With 
the exception of this energy interval and the interval 
from 30 to 55 MeV where statistics are very poor, the 
measured positron fraction lies between .24 and .35 from 
11 MeV up to ~300-400 MeV. Also shown in Figure 31 is 
the interstellar positron fraction for the collision 
source. We have derived this curve by combining the cal- 
culation of Ramaty and Lmgenfelter (28) for the pion-decay 
source with those of Abraham, Brunstexn, and Cline (30) 
for the knock-on source. In Table 6 we use the calculated 
and measured positron fractions to determine the relative 
contributions of directly accelerated (primary) and colli- 
sion-produced (secondary) electrons to the total cosmic-ray 
electron spectrum. We assume here and m the subsequent 
discussion that the collision source is the only signifi<- 

cant source of positrons of these energies. Energy loss 

1 

during the modulation process,, which we have ignored, 
affects the derived primary and secondary electron contri- 
butions. We shall return to ( this point later. 
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Fig. 31. Positron fraction of the extraterrestrial 

electron flux. N . / (N , + N ). The dashed 

e+ e+ e~ 

curve is the interstellar positron frac- 
tion for the collision source alone. It 
is derived by the author by combining the 
work of Ramaty and Lingenfelter (Ref. 28) 
and Abraham, Brunstein, and Cline 
(Ref. 30). 






TABLE 6 


POSITRON FRACTIONS 


Measured energy 
interval (MeV) 

Measured positron 
N e+ 

fraction F m = ^ — 
m N e + N e - 

Positron fraction for 
pion-decay source F^ (1) 

Relative contribution of 
pion-decay secondaries 
to measured electron 
flux m _ _B. (2) 

J, TT F 

TT 

Positron fraction for 

secondary electron 

source F (3) 
s 


11.1-17.1 17.1-30.2 

. 35± . 08 . 28± . 11 

.9 .9 

. 39± . 09 . 31± . 12 

.65 .9 


30.2-55.2 55.2-104.9 

. 12± . 24 

.9 .85 

. 13± . 27 

.9 .85 


Relative contribution 

of primaries to measured .46±.12 .69±.12 .87±.27 

electron flux 



( 2 ) 


(1) Ramaty and Lingenfelter (Ref. 28). 

(2) Derivation of the listed values assumes insignificant energy loss during 

(3) Includes both the interstellar pion-decay (Ref. 28) and knock-on (Ref. 3 


104.9-203.9 
. 24± . 14 

.85 

. 28± . 16 
. 85 

. 72+.16 


modulation, 
i) sources. 
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In Figure 32 we plot our measured positron flux 
together with values reported by Fanselow, Hartman, Hilde- 
brand, and Meyer (35) , Kniffen, Cline, and Fichtel (44) , 
and Cline and Porreca (43) . The solid curve, e + , is the 

5 

interstellar equilibrium positron spectrum from pion decay 
as calculated by Ramaty and Lingenfelter (28). We define 
the modulation factor for a given energy interval as the 
ratio of the measured intensity to the calculated inter- 
stellar intensity for the same energy interval. The 
modulation factors which we derive are shown in Figure 33. 
This graph suggests rapidly decreasing modulation below 
about 50 MV. For example, the modulation factor between 11 
and 17 MV is a factor of 50 greater than that between 30 
and 55 MV (i.e., the modulation is less at the lower 
energies ) . 

Fanselow et al . have suggested as a good fit to 
their own data the modulation function 

exp (-600/3R) for R > R = 300 MV 

° (4) 

exp (~600/3R o ) for R < R Q = 300 MV 

where the particle rigidity R is measured in MV. Use of 
eq. (4) with e + in Figure 32 yields curve a in Figures 32 
and 33. This functional formv is predicted by the 
diffusion-convection theory of solar modulation without 


3 E (R, 1965 ) 
3co(R) 



-138- 


Figure 32. Differential kinetic-energy spectrum 
of extraterrestrial positrons, inclu- 
ding data of other investigators. 

The year in which each measurement 
was made is indicated. The solid 
curve, e"t, is the interstellar posi- 
tron spectrum from pion decay ‘ 
calculated by Ramaty and Lingenf elter 
for an integral path length of 4 g/cm 2 
(Ref. 28). The modulated spectra 
a, b, and c are described m the 


text. 



POSITRONS /(M 2 -SEC-SR-MeV) 



KINETIC ENERGY (MeV) 


Figure 32 
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RIGIDITY (MV) 

Fig. 33. Absolute solar modulation factors, 
j E (R, t) /j^ (R) , for positrons in 

1968. Curves a and b are modulation 
functions described in the text and 
also shown in Figure 32. 
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energy loss (equations (1) and (3)). The parameters are 
xn reasonable agreement with values derxved from data for 
cosmic-ray deuterium and helium-3 (91) . On the other hand, 
curve a represents a very poor fit to our own data, having 
a chi-square probability much less than one percent. Due 
to the large error limits on our measurements, however, it 
is not possible to entirely exclude a similar modulation 
function with somewhat lower R , e.g., curve b m Figures 
32 and 33 which is 


j E (R, 1968) 
3co(R) 


exp (-600/PR) R > R q = 175 MV 

exp (-600/3R_) R < R = 175 MV 

u — o 


(5) 


and has a chi-square probability of 40 percent. This value 
of R 0 would also be in approximate 1 agreement with the 
rigidity at which Jokipii (82) suggests a change to 13- 
dependent modulation on the basis of the observed power 
spectrum of the interplanetary magnetic field. Note that 
eq. (5) is an equally good fit to the data of Fanselow, 
et al . In the following discussion we shall take the 
derived modulation factors at face value; the fact that 
modulation functions such as eq. (5) , which gives constant 
modulation at low rigidities, cannot be excluded should 
be kept in mind, however. 
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Our measured positron fractions (Figure 31) can be 
used together with the calculated interstellar positron 
spectrum (curve e* in Figure 32) to derive the galactic 
electron spectrum. In order to do this, however, it is 
necessary to know how much energy is lost by the electrons 
in reaching the inner solar system, i.e., at what energy 
m interstellar space the measured positron fractions are 
to apply. If we assume the energy loss to be an insigni- 
ficant fraction of the original energy, we obtain curve 2 
m Figure 34 for the galactic electron spectrum. The 
spectrum has been smoothly extrapolated above 200 MeV 
assuming no modulation above 5 GeV. Use of our measured 
positron fractions assuming no energy loss is equivalent 
to assuming equal modulation for the positron and negatron 
components. Also shown in Figure 34 are three estimates 
of the interstellar electron spectrum made by Webber (92), 
Verma (93), and Anand, Daniel, and Stephens (94) based on 
measurements of the galactic radio background. Although 
there is little actual overlap with our directly derived 
values, these estimates agree reasonably well with our 
extrapolated spectrum, particularly the estimate of 
Webber . 

Curve 3 in Figure 34 is the calculated interstellar 
secondary electron spectrum for 4 g/crn^ integral path 
length including both ,pion-decay (28) and knock-on 
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Fig. 34. Differential kinetic-energy spectra of inter- 
planetary and of interstellar electrons. 

Curve 1: composite interplanetary electron 

spectrum observed near the earth 
(see Figure 30) . The closed 
squares are our measurements. 

Curve 2: interstellar electron spectrum 

derived from curve 1 assuming 
insignificant energy loss during 
modulation. Our demodulated data 
points are shown as open squares. 

Curve 3 : interstellar secondary electron 

spectrum resulting from pion- 
decay (Ref. 28) and knock-on 
processes (Ref. 30). 


Curve 4: interstellar primary electron 

spectrum, obtained by subtracting 

i 

curve 3 from curve 2. 


Also shown are galactic electron spectra ob- 
tained by Webber (Ref. 92 ), Verma (Ref. 93), 
and Anand, Daniel, and Stephens (Ref. 94), 

I 

from studies of the galactic background syn- 

1 

chrotron emission. 



DIFFERENTIAL FLUX (ELECTRONS/M 2 - SEC -SR ■ MeV) 



KINETIC ENERGY (MeV) 


Figure 34 
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processes (30) . Subtraction of curve 3 from curve 2 yields 
curve 4 for the equilibrium spectrum of primary electrons 
in the galaxy. This spectrum appears to have a broad 
peak aroung 30-50 MeV which, if real, might be associated 
with a minimum escape energy from the source region. 

Curve 1 in Figure 34 is a composite spectral shape for 
electrons measured near the earth in 1967-1968. Below 
about 10 MeV this spectrum is primarily the work of 
Simnett and McDonald (21) . These authors have noted pre- 
viously that their measured spectrum is compatible with a 
sole origin in the galactic knock-on component (the 
dominant contributor to the collision source below about 
20 MeV) provided solar modulation is insignificant at 
these energies. Recently Beedle, Lezniak, Rockstroh, and 
Webber (95) have reported electron measurements which indi- 
cate that the similarity to the galactic knock-on spectrum 
may persist down to 200 keV. Our data imply that the 
modulation, though high near 100 MeV, does, in fact, seem 
to be rapidly decreasing below about 50 MeV. Furthermore, 

in Figure 32, extrapolation of the calculated positron 

( 

spectrum, e*, to lower energies would intersect the two 
lowest values of Cline and Porreca, which are thus consis- 
tent .with no modulation. Their value between 5 and 10 MeV 
is inconsistent with no modulation, however. We note, 
however, that the spectral shape reported by Cline and 
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Porreca is one of several interpretations of their data. A 
mean flux value for the entire interval from 2 to 10 MeV 
is perhaps equally valid; such a value would fall reasona- 
bly close to the interstellar positron intensity from pion 
decay (96) . The absence of significant modulation at low 
energies would also help to explain the lack of long-term 
variations in the 3-12 MeV election flux which has been 
noted by Cline and McDonald (20) . This absence is espe- 
cially striking m view of large short-term variations of 
this component seen by the same authors. The interplane- 
tary magnetic field has been observed to have a filamen- 
tary structure (97-99) with flux tubes down to perhaps 
4000 km m diameter (approximately the cyclotron radius of 
a 70 MV particle). Beuermann, Rice, Stone, and Vogt (100) 
have noted that a well ordered field within such flux 
tubes might enable low-rigidity particles to penetrate 
into the inner solar system without undergoing significant 
modulation. Parker (83) has also suggested, on the basis 
of measurements of the interplanetary magnetic field, that 
low-energy electrons (4 10 MeV) may have relatively free 
access to the inner solar system. 

Up to this point we have not considered the adia- 
batic cooling of the cosmic-ray particles diffusing through 

I 

the expanding solar wind. The energy loss due to this 
effect may not be insignificant, however. Parker (84) has 


I 
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estixnated that a cosmic-ray particle which penetrates to 
the earth will have lost 5-20 percent of its original 
energy. At the low rigidities with which we are primarily 
concerned the energy loss could be even higher, however, 
perhaps as high as 50 percent or more. We shall examine 
here the general effect which energy losses of this magni- 
tude would have on our earlier discussion. 

The relative contribution, T^, of pion-decay secon- 
daries to the electron flux measured at energy E is 


T {E, E 1 ) 
ir 


(E’) 


where F M is the positron fraction measured at the earth, 

(EM is the positron fraction calculated for the pion- 

decay source at energy E' in interstellar space, and 

AE = E' - E is the energy lost by the electrons in reaching 

the earth. In Table 6 we list the values of T derived 

7T 

for each of our 5 energy intervals assuming AE = 0. It is 

apparent from Table 6 that both and F^ are relatively 

constant from 11 to 200 MeV. It follows that the derived 

values of do not significantly change for AE/E' as high 

as .5 or qreater. This is true whether T is considered 

rr 

at the earth (i.e., at fixed E) or in interstellar space 
(i.e., at fixed EM* A direct consequence is that phe 
derived interstellar, differential kinetic-energy spectra 
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of all electrons and of directly accelerated electrons, 
curves 2 and 4 in Figure 34, are not appreciably altered 
even for quite a large energy loss (AE/E' as large as 
.5-. 7). Because of adiabatic deceleration, the derived 
interstellar spectra will not, of course, reach to energies 
as low as those we have measured at the earth. Furthermore 
our lowest measured interval, 11.1-17.1 MeV, contains an 
insignificant knock-on contribution for AE/E' as small as 
.2-. 3'. Thus, even though there would still be a good 
indication of a considerable flattening in curves 2 and 4 
in Figure 34 below several hundred MeV, support for the 
sole knock-on origin of electrons below ~10 MeV must be 
considered weakened. 

Adiabatic cooling of the cosmic-ray particles may 
itself cause some of the apparent decrease in solar modu- 
lation which we observe at low positron energies. Since 
no particles are actually lost due to adiabatic decelera- 
tion, the presence of the broad maximum around 30-35 MeV 
m the interstellar differential kinetic-energy spectrum 
of positrons (e^ in Figure 32) , coupled with a significant 
energy loss, leads to an accumulation of positrons at low 
energies. We shall illustrate the approximate effect by 

I 

including this energy loss with the simple modulation 
function given by eq. (5) . The discussion will be similar 
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to a section in Parker (83) with some change of notation. 
Derivations of the relevant equations can be found in this 
source or references quoted therein. 

We shall consider a radial solar wind with constant 
velocity v and let the cosmic-ray diffusion coefficient «■ 
be isotropic and uniform out to the boundary of the modu- 
lation region at r - L. In the absence of energy loss, 
the cosmic-ray intensity at the earth, j E (R,t), is reduced 
relative to the interstellar intensity j^CR) by the factor 

3e (R't) 

bluo — = exp (_Lv/K) <7) 

Comparison of eqs . (5) and (7) shows that the modulation 

function assumed m eq. (5) is equivalent to 

R > R q = 175 MV 

( 8 ) 

R <_ R Q = 175 MV 

The rate of energy loss of a relativistic particle 
which is undergoing adiabatic deceleration is given by 

1 dE 2v 1 

E dt 3r t E 

where t E is the characteristic time of the energy loss 


Lv 

600 

[ 600 PR o 


rate. 
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We now make two further approximations. First, we 

assume that as a particle loses energy, the diffusion 

coefficient which applies to its motion does not change, 

i.e., k is constant for a given particle and, in particular, 

is the value applicable to the rigidity R' of the particle 

when it first entered the modulation region. With the 

assumed form for k (eq. (8)), this condition is fulfilled 

in any case for R'^175 MV. Thus the low-rigidity interval 

with which we are primarily interested is not affected. 

Second, we substitute for r in eq. (9) a characteristic 

value r such that t_ can be taken to be a constant. With 
o E 

these s implications , it follows directly from integrating 
eq. (9) that the energy of a cosmic-ray particle at the 
earth is 


E = E' exp (- t L /t E ) (10) 

Here t L is the time required to diffuse into the inner 
solar system. It is approximately 

t_ = L 2 /4K (11) 

-L 

Combining eqs. (8), (10), and (11), we obtain 
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E 



L E ' 



exp 



100 L 


3R'r, 


100 L 


o J 




R' > R = 175 MV 
o 

R’ < R = 175 MV 
— o 


(12) 


We can obtain an estimate of the energy loss of the cosmic- 
ray electrons by substituting r Q = L into eq. (12). In 
that case we find, for example, that E - . 6E’ for E'£l75 
MeV and E ^ . 9E 1 for E' = 1000 MeV. 

The total effect on the interstellar positron 

spectrum e* is derived by including energy loss as given 
□ 

by eq. (12) with the simple diffusion-convection modulation 
as given by eq. (5) . The energy loss is normalized by 
taking E = . 5E' for E'<175 MeV and we substitute R’ for R 
on the right-hand side of eq. (5) (equivalent to an 
unchanging diffusion constant k for a given particle) . In 
this manner we obtain curve £ in Figure 32 for the differ- 
ential kinetic-energy spectrum of positrons at the earth. 
Curve c is to be compared with curve b which follows from 
eq. (5) without energy loss. The accumulation of particles 
at low energies is apparent; the greater similarity of 
curve £ to the shape of our measured spectrum is also evi- 
dent. 

A measurement of the positron spectrum by itself 
does not lead to a determination of the magnitude of the 
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adiabatxc deceleration. It seems likely that adiabatic 
deceleration has played some role in shaping the inter- 
planetary positron spectrum which we observe. On the 
other hand, the general arguments advanced m support of 
decreasing modulation at low rigidities seem to result 
an a self-consistent picture which correlates many differ- 
ent observations. It is possible that both processes, 
adiabatic deceleration and decreasing modulation at low 
rigidities, contribute to shaping the cosmic-ray electron 
spectrum which we observe at the earth. 



VI. CONCLUSIONS 


In this thesis we have presented measurements of the 
cosmic-ray positron and negatron spectra between 11 and 204 
MeV. The data lead to several conclusions. 

1) We obtained new data relating to the behavior of 
charged particles inside the earth's magnetosphere. Our 
measurement of the diurnal intensity variation of electrons 
indicates that the local geo-magnetic cutoff at invariant 
latitude A = 68°-70° was <11 MeV at night. The approxi- 
mately equal contributions of positrons and negatrons in 

the daytime electron flux (Figure 27) is in agreement with 
the charge ratio expected for return albedo particles and 

is in marked contrast to the measured nighttime charge 

ratio. This feature supports the interpretation of the day- 
time electrons as return albedo and also allows us to eli- 
minate albedo as a significant contamination in the night- 
time electron flux. 

2) We derived interplanetary positron and negatron 
spectra (Figure 29) . Our total electron (e + +e~) spectrum 
agrees well with recently reported satellite measurements 
(Figure 30) . We measured a positron- to-total-electron ratio 
which lies between .24 and .35 over the major part of our 
energy interval. By comparing our measured positron flux 
with a calculated interstellar positron spectrum resulting 
from pion decay (28) , we derived the magnitude of the abso- 
lute solar modulation of positrons in 1968 (Figure 33) . 

The measured positron intensity near 50 MV was less than 
~.03 times the interstellar flux; near 15 MV, however, the 
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modulation factor was -.25. This apparent absence of sig- 
nificant modulation at low energies supports the suggestion 
that the knock-on component contributes significantly to 
the interstellar electron flux below -10 MeV (21) . 

3) We derived the total interstellar electron spec- 
trum, ignoring adiabatic deceleration during modulation 
(Figure 34) . The extrapolation of this spectrum to higher 
energies agrees reasonably well with estimates based on the 
galactic background radio emission. We also derived the 
interstellar spectrum of directly accelerated (primary) 
electrons. This spectrum exhibits considerable flattening 
below several hundred MeV when compared with the electron 
spectrum at higher energies . We examined the general 
effect of adiabatic deceleration during the modulation 
process and found that many of our deductions (in particular, 
the derived interstellar total electron spectrum and the 
directly accelerated component) are not significantly 
affected, even for an energy loss as high as 50 percent 
of the original energy. 
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